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PREFACE 


I gave the lecture series "Boundary— Layer Theory" in the winter 
semester 1941/42 for the members of my Institute and for a considerable 
number of collaborators from the Hermann Goring Institute for Aviation 
Research. The series embraced a total of sixteen two-hour lectures. 

The aim of the lecture series was to give ,a survey of the more 
recent results of the theory of viscous fluids as far as they are of 
importance for actual applications. Naturally the theory of the boundary 
of frictional layer takes up the greatest part. In view of the great 
volume of .material, a complete treatment was out of the question. 

However, I took pains to make concepts everywhere stand out clearly. 
Moreover, several important typical examples were treated in detail. 

Dr. H. Hahnemann (LFA, Institute for Motor Research) went to 
considerable trouble in order to perfect an elaboration of this lecture 
series which I examimed and supplemented in a few places. Miss Hildegard 
Miinz participated in the illustration. To both I owe my most sincere 
thanks for this collaboration. 


Schlichting 


Aerodynaml sches_ Institut 

der Technischen Hochschule, Braunschweig 

October 1942. 
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TECHNICAL MEMORANDUM NO- 1217 


LECTURE SERIES "BOUNDARY LAYER THEORY" 
Part I — Laminar Flows* 

By H. Schlichting 


First lecture (Dec. 1, 1941) 


INTRODUCTION 


Gentlemen: In the lecture series starting today I want to give 

you a survey of a field of aerodynamics which has for a number of years 
been attracting an ever growing interest. The subject is the theory of 
flows with friction, and, within that field, particularly the theory of 
friction layers, or boundary layers. 

As you know, a great many considerations of aerodynamics are based 
on the so-called ideal fluid, that is, the frictionless incompressible 
fluid. By neglect of compressibility and friction the extensive mathe- 
matical theory of the ideal fluid (potential theory) has been made 
possible. 

Actual liquids and gases satisfy the condition of incompressibility 
rather well if the velocities are not extremely high or, more accurately, 
if they are small in comparison with sonic velocity. For air, for 
instance, the change in volume due to compressibility amounts to about 
1 percent for a velocity of 60 meters per second. 

The hypothesis of absence of friction is not satisfied by any 
actual fluid; however, it is true that most technically important fluids, 
for instance air and water, have a very small friction coefficient and 
therefore behave in many cases almost like the ideal frictionless fluid. 

Many flow phenomena, in particular most cases of lift, can be treated 
satisfactorily, — that is, the calculations are in good agreement with 
the test results, — under the assumption of frictionless fluid. However, 
the calculations with frictionless flow show a very serious deficiency; 
namely, the fact, known as d’Alembert’s paradox, that in frictionless 
flow each body has zero drag whereas in actual flow each body experiences 
a drag of greater or smaller magnitude. For a long time the theory has 
been unable to bridge this gap between the theory of frictionless flow 
and the experimental findings about actual flow. The cause of this 
fundamental discrepancy is the viscosity which is neglected in the theory 

♦"Vortragsreihe **Grenzschichttheorie. ’ Teil A: Laminare Stromungen." 

Zentrale fur wissenschaftliches Berichtswesen der Luftfahrtforshung des 
Generalluftzeugmeisters (ZWB) Berliit-Adlershof , pp. 1—153* Given in the 
Winter Semester 1941/42 at the Luftfahrtforschungsanstalt Hermann Uoring, 
Braunschweig. The original language version of this report is divided into 
two mw.in parts, Teil A and Teil B, which have been translated as separate 
NACA Technical Memorandums, Nos. 1217 and 1218, designated part I and part II, 
respectively* 
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of the ideal fluid; however* in spite of its extraordinary smallness it 
is decisive for the course of the flow phenomenon. As a matter of fact 
the problem of drag can not he treated at all without taking the viscosity 
into account. 

Although this fact had been known for a long time* no proper approach 
to the theoretical treatment of the drag problem could be found until the 
beginning of the present century. The main reason was that unsurmountable 
mathematical difficulties stood in the way of theoretical treatment of the 
flow phenomena of the viscous fluid. It is Professor Prandtl' s great merit 
to have shown a way to numerical treatment of viscosity* particularly of 
the technically important flows under consideration and thereby to have 
opened up new viBtas on many important perceptions about the drag problem 
and related questions. Prandtl was able to show that in the case of most 
of the technically important flows one may treat the flow* as a whole* as 
frictionless and utilize the simplifications for the calculation thus made 
possible* but that in the immediate neighborhood of the solid walls one 
always had to take the friction into consideration. Thus Prandtl subdivides, 
for the purpose of calculation* the flow surrounding a body into two 
domains: a layer subject to friction in the neighborhood of the body* and a 

frictionless region outside of this layer. The theory of this so-called 
”Prandtl's friction or boundary layer” has proved to be very fruitful in 
modern flow theory; the present lecture will center around it. 

At this point I want to indicate a few applications of the boundary- 
layer theory. A first important application is the calculation of the 
frictional surface drag of bodies immersed in a flow* for instance, the 
drag of a flat plate in longitudinal flow* the frictional drag of a ship* 
a wing profile, and an airplane fuselage. A special property of the 
boundary layer is the fact that under certain circumstances reverse flow 
occurs in the immediate proximity of the surface. Then* in connection 
with this reverse flow* a separation of the boundary layer takes place* 
together with a more or less strong formation of vortices in the flow 
behind the body. Thus a considerable change in pressure distribution, 
compared with frictionless flow* results* which gives rise to the form 
drag of the body immersed in the flow; The boundary— layer theory therefore 
offers an approach to the calculation of this form drag. Separation occurs 
not only in the flow around a body but also in the flow through a divergent 
tunnel. 

Thus flow phenomena in a diffuser, as* for instance* in the bucket 
grid of a turbine, may be included in boundary— layer theory. Furthermore* 
the phenomena connected with the maximum lift of a wing* where flow 
separation is concerned* can be understood only with the aid of boundary- 
layer theory. The problems of heat transfer also can be explained only 
by boundary— layer theory. 

As will be shown in detail later* one must distinguish between the 
two states of boundary— layer flow — laminar and turbulent; their flow laws 
are very different. Accordingly* the lecture will be divided into three 
main parts: 1. Laminar flows* 2. Turbulent flows, 3* Laminar— turbulent 

transition. Although the boundary layer will be our main consideration* 
it will still be necessary as preparation to discuss to some extent the 
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general theory of the viscous fluid. This will he done in the first 
chapter. 


CHAPTER I. VISCOSITY 


Every fluid offers a resistance to a form variation taking place in 
finite time interval, which is of different magnitude according to the 
type of fluid. It is, for instance, very large for syrup or oil, but 
only small for the technically important fluids (water, air). 

The concept of viscosity can be best made clear by means of a test 
according to figure 1: 

Let fluid be between two parallel plates lying at a distance h 
from each other. Let the lower plate be fixed, while the upper plate is 
moved with the velocity Uq uniformly and parallel to the lower one. 

For moving the upper plate a tangential force P must be expended which 
is 


p = uF — 
h 


( 1 . 1 ) 


according to experiment, where F Is the area of the upper plate and 4 
is a constant of proportionality. (End effects are not included). The 
quantity p is called the viscosity coefficient or the dynamic viscosity. 

Since the phenomenon in question is a parallel gliding, the transverse 
velocity component in the y-direction, denoted by v, equals zero. The 
fluid adheres to the upper and lower surface, respectively, a linear 
velocity distribution between the plates is set up, the magnitude of 
which depends solely on y. 


/ v y du u o 

u(y) = “o S or: £ = T 


Since for y = 0: u = 0, for y = h: u = Uq. If one designates P/F, 

the tangential force per unit area, as the frictional shearing stress 
there follows: 


= 4 7^ kg/m2 


T 


( 1 . 2 ) 
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Ths dimensions of n are accordingly kg sec/mt A flow as represented 
in figure 1, where no transverse velocity occurs and the shearing stress 
at all points of the flow is therefore given by equation -(1.2), is called 
simple shear flow. In the special case described, the shearing stress is 
everywhere of equal magnitude, and equal to that at the surface. Besides 
the dynamic viscosity |i the concept, of kinematic viscosity v is 

required, which for the density pU sec 2nr-4~] 

is defined as 


V = H [m a /a] 


For 20° C, v is, for instance, for water: 


V 


= 1.01 X 10 


,-6 



for air: 


v 


= 14.9 x 10 -6 m 2 /s ~ 


- 10”^ m 2 /s 

7 / 


if the air pressure has the standard value p Q = 760 mm hg. 

CHAPTER II. POISEU ILLE FLOW THROUGH A PIPE 


The elementary empirical friction law of the simple shear flow 
derived above permits the immediate determination of the flow and the 
resistance in a smooth pipe of circular cross section and of constant 
diameter, d = 2r. At a very large distance from the beginning of the 
pipe one cuts off a piece of pipe of length l (fig* 2) and examines 
the cylinder of diameter 2y, the axis of which is identical with the 
pipe axis. According to what has been said so far, the velocity probably 
will be again a function of y. A pressure difference p^ — p£ is 

required for forcing the fluid through the cylinder. According to 
practical experience, the static pressure across every cross section may 
be regarded as constant. The flow is assumed to be steady and not 
dependent on the distance from the beginning of the pipe. Equilibrium 
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must then exist between the pressure and. the frictional shearing stress 
which attempts to retard the motion. Thus for the cylinder of radius y 
the following equation is valid: pressure force difference acting at 

the cross sections = frictional force acting along the cylinder wall, or 



Jty2 = 2itylT 


( 2 . 1 ) 


or 


Pi ~ P2 y 

l 2 


(2.1a) 


Sine© flow parallel to the axis is to be expected, one takes from the 

previous paragraph, t = - |i— (the minus sign indicates that the 

dy 

velocity diminishes with increasing distance from the axis; thus du/dy 
is negative, the shearing stresses under consideration, however, are 
positive), and, after separation of the variables, du becomes: 


du = — 


Pi ~P 2 

ixl 



( 2 . 2 ) 


and, on integration: 


u(y) 


Pi ~ P2 




From the fact that for y = r the velocity is supposed to be u(y) = 0 
follows that the constant of integration C has to be C = r2 /k . Thus : 


- Sfcr 2 K - * 2 ) 


(2.3) 


6 


NACA TM No. 1217 


This equation (2.3) is Poiseuille’s law for pipe flow. It states that 
the velocity u(y) is distributed parabolically over the pipe cross 
section. The apex of the parabola lies on the pipe axis; here the 
velocity is greatest, namely: 


u max 


h ~ p 2 rP 


Therewith one may write (2.3): 


= w ( x - 4) 


(2.4) 


(2.3a) 


By Poiseuille’s law (equation (2.3)) the drag of the developed laminar 
flow (which is proportional to p 1 — p 2 ) is directly proportional to 

the first power of the velocity. 

This statement is characteristic of all kinds of laminar flow 
whereas, as will he seen later, the drag in turbulent flow is almost 
proportional to the second power of the velocity. 


The flow volume for the present case remains to he given. With 

n 


designating an area element, Q 
velocity paraboloid, therefore 


is Q = 


u(y) d F = volume of the 


0 


dF 


Q = i nr 2 u 

2 max 


jtr 


*1 (P 1 ‘ 


(2.5) 


This flow law is often used for determination of the viscosity, by 
measuring the quantity flowing subjected to a pressure gradient (usually 
produced by gravity in a vertical capillary tube). Of course, the 
starting losses must be taken into consideration which due to the mixing 
zone (vortex formation) at the pipe end are not recovered to their full 
extent. 

A drag coefficient X will now be defined. Since turbulent flows 
are more important than laminar ones and since the drag in turbulent 
flow increases about as the square of the velocity, X will also be 
referred to u^. 

For flow problems, let X thus be defined as: ratio of the pressure 

drop along a test section of a specified characteristic length to the dynamic 
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pressure q = pu^/2, with u = - 0 — = the mean velocity (average taken 

nr 2 

across the cross section). Then: 


§£. _A_ 

dx rr 


( 2 . 6 ) 


with d = characteristic length, thus, for the present case, the pipe 
diameter, and with \ = dimensionless quantity. For the present developed 
laminar pipe flow, according to equation (2.5) 


dp = Pi ~ p 2 = 8uQ 


Thus: 


= 2d_ 8pQ = 2d8n Q = l6dp = 64u/o 


_2 4 

pu irr 


— 2 2 _ 

pur nr u pufj^j 


(2.7) 


with the dimensionless quantity Re = — signifying the Reynolds number 

of the circular pipe. Since the pressure drop which is only linear ily 
dependent on the velocity was referred to u^, then, for laminar flow; 

A ~ A logarithmic plot of A = f(Re) or A = f(u) therefore results 
in a straight line inclined 45° toward the Re— axis (compare fig. 82 Part II.) 


After this short analysis of the one-dimensional case of viscous 
fluid we will now consider the three-dimensional case. 
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CHAPTER III. EQUATIONS OF MOTION OF THE VISCOUS FLUID 

a. State of Stress 


For this purpose one must know first of all the general state of 
stress in a moving viscous fluid and must then connect this state of 
stress with the state of deformation. For the deformation of solid 
todies the resistance to the deformation is put proportional to the 
magnitude of the deformation (assuming the validity of Hooke's law). 

For flowing fluids, on the other hand, the resistance to deformation 
will depend on the deformation velocity, that is, on the variation of 
velocity in the neighborhood of the point under consideration. (Solid 
bodies: displacement gradient = displacement per second. Fluid: 

velocity gradient). 

One starts from the basic law of mechanics according to which: 
mass x acceleration = sum of the acting, or resultant force. For the 
mass-per— unit volume, that is, the density p, one may write the law 


p *D* 

Dt 


= K + R + F 


(3-1) 


Dw 

— = substantial acceleration 
Dt 

K = mass forces 

R = surface forces, composed of pressure forces normal to the 
~ surface and frictional forces in the direction of the 

surface 

F = negligible extraneous forces 


In order to formulate the surface forces, one imagines a small rectangular 
element of volume dV = dx dy dz cut out of the flow (fig. 3) the left 
front corner of which lies at the point (x, y, z). The element is to be 
very small so that only the linear variations of a Taylor development 
need to be taken into consideration; on its surfaces dy dz act the 
resultant stresses (vectors): 


J?x or £x 



respectively 


(3.2) 


* Throughout the text, underscored letters are used in place of 
corresponding German script letters used in the original text. 
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(The index x signifies that the stress tensor acts on a surface element 
normal to the x— direction) . 


Analogous terms result for the surfaces dz dx normal to the 
y-axis and dx dy normal to the z-axis, if x in equation (3.2) is 
replaced everywhere hy y or z, respectively. From this there results 
as components of the resultant force: 

dp x 

Force on the surface element normal to the x— direction: dx dy dz 

3x 


dp. 


Force on the surface element normal to the y— direction: — — dy dz dx 

dy 


dp 


Force on the surface element normal to the z-direction: — - dz dx dy 

dz 


The total resultant surface force E per unit volume caused hy the 
state of stress is therefore: 


E 


dx dy dz 


(3-3) 


p , p and p are vectors which can he further decomposed into 

—x — y — z 

components. In this decomposition the components normal to every surface 
element, that is, the normal stresses, are designated hy a (indicating 
hy the index the direction of this normal stress); the other components 
(tangential stresses) are denoted hy t (with double iqdex: the first 

indicates to which axis the surface element is perpendicular, the second, 
the axial direction of the stress t). With these symbols there is: 


p = ic + Jr + kT 
-x x xy xz 


< p = iT + jg + kT > 

t-y- yx y yzf 


p = It + Jt + kcr ! 
- z zx zy z 
^ J 


(3.4) 


This state of stress represents a tensor with nine vector components, 
which can he characterized hy the stress matrix (stress tensor): 
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xy 


zy 



(3-5) 


It can readily "be shown that those of the six tangential stresses 
which hare the same indices, although in interchanged sequence, must be 
equal. This follows for a homogeneous state of stress from the equilibrium 
of the small cube dx dy dz with respect to rotation: 

Since r dy dx is the force attempting to rotate the cube counter— 
xy 

clockwise about the z-axis, (seen from above in fig. 3 ) > with the lever 
arm dx, and since, correspondingly, the force — t dx dz attempts to 

rotate the cube clockwise about the z-axis, with the lever arm dy the 
balance of moments requires: 


dy dz dx - T xy dx dz dy = 0, thus T xy ■ t^. 


Correspondingly, because of freedom from rotation about the x— axis 

r = T , and because of freedom from rotation about the y-axis 
yz zy* 

t = t , the nine components of the stress tensor are reduced to six 
zx xz* 

and the stress matrix (equation ( 3 . 5 )) is converted into the stress 
matrix symmetrical with respect to the principal diagonal: 



(stress matrix) (3*6) 


For the frictional force one obtains according to equation (3*3) by 
insertion of the components from equation ( 3 -h) and by reduction to the 
six remaining terms according to equation ( 3 . 6 ): 
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R = i 


St St \ ( St So St 

x xy ' xz . f x.j y yz 
+ + 1 + J I + — - + — - 

Sx Sy Sz/ \ Sx Sy Sz 


= x— component 


+ y— component 


St. 


St„ 


+ k 


xz + ~ ’yz + 
Sx Sy 


Sa 2 

Sz. 


(3.7) 


+ z— component. 

For the case of the frictionless (ideal) fluid all shearing 
stresses disappear 


r = r = r = 0 
xy yz zx 


(3.8) 


and only the normal stresses remain, which in this case are all equal. 
Since the normal stresses from within toward the outside are denoted as 
positive, the normal stresses equal the negative fluid pressure: 

a x = a y = a z = - p (3>9) 

The static pressure equals the negative arithmetic mean of the normal 
stresses: 

_p = i/(j + a + a \ 

* 3 \ x y V 


b. State of Deformation 

The state of stress treated so far is, alone, not very useful. 
Therefore we will now consider the state of deformation (that is the 
field of velocity variations) and then set up the relations between 
state of stress and state of deformation. 
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Let the velocity with the components u^, v A , in the 

directions of the axes exist at the point A the coordinates of which 
are x A , y A , z A> 


If one limits oneself to the points x, y, z in th9 immediate 
neighborhood of A with the velocity w = iu + jv + kw, and if one 
limits oneself — as also in setting up the state of deformation — to 
linear terms only, one obtains for the deformation the relative change 


in position between the points x, y, z and x 


A* 




z per unit 


time, that is, the difference of the velocities at the points x, y, z 



J 


dw = i du + J dv + k dw = distortion of the fluid region in the 
neighborhood of the point A. 


Omitting the index A one obtains therefore: 


, dw „ dw„ <3w „ . / du . du , du . \ 

dw = — dx + -=-dy + -^dz = i( — dx + — dy + — dz) 

ax ay oz \ax ay dz / 

(3.11a) 

+ J (^ ai + ^ dr + sr d 0 + £ fe to + ¥ d:r + af 4z ) 


Thus the velocity variation (and hence, on integration, the velocity 
itself) in the neighborhood of the point A is known if the nine 
partial derivatives of the velocity components with respect to the space 
coordinates are known. Corresponding to the stress matrix, one may form 
a deformation matrix: 
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du 

<3u 

du 

Six 

dy 

dz 

dv 

dv 

dv 

dx 

dy 

dz 

dvr 

d w 

dw 

dx 

dy 

dz 



deformation matrix 


(3*12) 


The friction forces of the viscous fluid are given "by a relation (which 
will have to be determined) between these two matrices. First, the 
deformation matrix is to be somewhat clarified. 

1. Case of pure elongation. 


One assumes u - u. = a(x — x A ), with a = ^ = constant. Let all 

“ A dx 

other terms of the matrix disappear; the matrix will then appear as 
follows : 



0 

0 

0 


0 

0 

0 


\ 

/ 


Then the velocity variation is simply du = a dx, and u = ax. 

All points of the y— axis remain at rest, the points to the right and 
left of it are elongated or compressed, according to whether a > 0 
or a <0 (fig. 5). The equation u = ax therefore represents an 
elongation or expansion parallel to the x— axis. Corresponding relations 
apply for the other terms of the principal diagonal of the matrix. 

2. Case of pure translation 


All terms disappear; the matrix then reads: 


0 0 0 
0 0 0 


0 


0 


0 



Ik 
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In this case, which, as a matter of fact, should have "been mentioned 
first, u — u a = 0, du = 0; u = constant. The velocity component 

parallel to the x— axis is uniform (correspondingly for the other axes). 
3. Case of angular deformation. 


One assumes u — u = e(y — y ), with e = — = constant. All 
a 8- 

other terms equal zero, and the matrix reads: 


0 e 0 
0 0 0 
0 0 0 


du = e dy and u = ey; that is, all points of the x-axis retain their 
position; all points of the y-axis shift to the right (left), when e > 0 
(e < 0); for e > 0 the y— axis is rotated clockwise "by the angle c 
(because of the linearity). The y-axis is simultaneously elongated. The 
phenomenon in question is therefore a shearing (fig. 6), with tan e = e. 

Correspondingly there results for v — v = f (x — x.) and 

Q, A 



M = 



0 

0 

0 


0 

0 

0 


dv = f dx; v = xx 


All points of the y— axis retain their position; the points of the x— axis 
are rotated by the angle 6; tan 5 = f (fig. 7). Terms outside of the 
principal diagonal of the matrix result therefore in a deformation of 
the right angle with axis— elongation (shearing). The right angle between 
the x— and y— axes is, therefore, for e > 0 and f > 0 deformed by 


_ ov + ou _ y _ Deformation about the z— axis 
dx dy xy 


€ + 6 
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Correspondingly : 


zx 


— + — = deformation about the y— axis 
dz cbc 


7 = — + — = deformation about the x— axis 

yz by bz 


(The deformation angles are herein regarded as small so that the tangent 
may be replaced by the argument). 


c. Navi er— Stoke s Formulation for the Stress Tensor 


One now proceeds to relate the stress matrix (equation (3*6)) with 
the deformation matrix (equation (3*12)). The former is symmetrical with 
respect to the principal diagonal, but not the latter. However, one 
obtains a symmetrical deformation matrix by adding to equation (3. 12) 
its reflection in the principal diagonal. Furthermore, one first splits 
off the pressure p (contribution of the ideal fluid) from the stress 
matrix and sets the remaining stress matrix, according to Stokes, 
proportional to the deformation matrix made symmetrical: 


CT x T xy T xz 


p 0 0 

r a t 



xy y yz 


0 p 0 

t 

j 

T xz T yz °z 


0 0 p 



du 

du 

\bu 

£v 

bw 1 

dx 

by 

bz 


bx 

bx 

dx ■ 

dv 

bv 


+ [i 

du 


dw 

ctac 

by 

bz 


by 

by 

by 

dw 

dw 

dw 




bw 

dx 

by 

bz 


bz 

bz 

bz 


From equation (3.13) each stress component may be given immediately by 
coordinating the homologous parts of the matrices to each other. For 
instance : 


= — p+2|a — = static pressure + pressure due to 
dx velocity variation, or: 
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°x * - * * "y " - P + ^ °2 - - p + 

T = T = 

*3 7* \$x Zy) 

T — T = u/^Si 4. 

V zy \zj % z ) 

r = r = 

T xz zx ^ 

[ Thus for one-dimensional flow t™. = I 

2* dy 

Furthermore , there follows from equation ( 3 .I 3 ): 

3 (° I+ °y + M = - p + f‘‘(^ + | + a?) 

or 

3 (°i + °j * °z) * " p (3-14) 

“because 

|u + |l + |w = cL lTW = 0 
ox dy dz 


for the incompressible flows free of sources and sinks under considera- 
tion.* Thus for the viscous incompressible flow, as for the ideal fluid 
the pressure equals the arithmetic mean of the normal stresses. 

With these results the components of the friction force may be 
expressed according to equation ( 3 . 7 ) as follows: 


* The compressibility manifests itself as normal stress, since it can 
be interpreted as a pressure disturbance, for instance due to variation 
in density, which attempts to spread in all directions — considered 
infinite simally . 
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iz * 3l 3 * U 8 3y 2 3z 2 


+ n ( + d 2y + 

y^.2 dxdy dxdz/ 


or, since div w = 0, 


dp 

Sac 


R x = - ^ + MAu 


Ej = .| +1JiT k, 


E = — + iiAw 

z 3z 


in 


which Au = 2-a + & + & 


bx 2 by 2 bz 2 


(3-15) 


If one finally designates the mass forces hy K = p(iX + JY + kZ), and 
assumes the decomposition of the substantial derivative into a local 
and convective part as known from Euler * s equation, one obtains for the 
components of the equation of motion of the non— stationary, incompressible, 
and viscous fluid from equation (3*1): 


fdu + u dji + T du + w £u\ 

\dt dx by dz / 

f bv bv bv bv 

\ dt bx by bz 

(&■+ 

\ dt dx 




dp 

= pX + |jAu 

bx. 


pY - + pAv } 

o y 


+ V dw + w ^ 
by bz 


= pZ — + pAw 

OZ 


(3.16) 


In addition, the continuity equation 
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is used. Written in vector form, the Navier-Stokes differential equation 
and the equation of continuity read 


Dw 

Dt 


K — - grad p + VAw 


(3-18) 


div w = 0 


(3-19) 


Due to the friction terms, therefore, terms of the second order enter 
the differential equation. 

Boundary conditions are attached to these equations. If all friction 
terms on the right side are cancelled, that is V = 0, the differential 
equations Become equations of the first order and one Boundary condition 
is sufficient, namely the Boundary condition of the potential flow: 
v n = 0 on the Bounding walls. 

This means that the normal component v n of the velocity at the 

Bounding surface must disappear on the surface itself whereas the fluid 
still can glide parallel to the Boundary (tangential velocity v^ 

parallel to the surface ^ 0). 

For viscous flow where the differential equation is of the second 
order, two Boundary conditions are required, namely: 


v n = 0 and v^ = 0 (condition of no slip) (3-20) 


that is, the fluid must in addition adhere to the surface. 

Second lecture (Dec. 8, 19^1) 


CHAPTER IV. GENERAL PROPERTIES OF NAVIER— STOKES EQUATIONS 


These Navier— Stokes differential equations represent together with 
the equation of continuity a system of four equations for the four 
unknown quantities u, v, w, p. On the left side of the Navier— Stokes 
differential equations are the inertia terms, on the right side the mass 
forces, the pressure forces, and the friction forces. 

Since Stokes’ formulation is, of course, at first purely arbitrary, 
it is not a priori certain whether the Navier-Stokes differential 
equations describe the motion of a fluid correctly. They therefore 
require verification, which is possible only By way of experimentation. 
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Unfortunately, due to unsurmountable mathematical difficulties, a 
general solution of the differential equation is not yet known, that is, 
a solution where inertia and friction terms in the entire flow region 
are of the same order of magnitude. However, known special solutions 
(for instance, the pipe flow with predominant viscosity or cases with 
large inertia effect) agree so well with the experimental findings, 
that the general validity of Navier-Stokes differential equation hardly 
seems questionable. 

The plane problem: 

By far the greatest part of the application of Navier-Stokes 
differential equations concern "plane” cases, that is, the cases where 
no fluid flows in one direction. The velocity vector w is then given 


w = iu (x, y, t) + Jv (x, y, t) 


(4.1) 


since w = 0. The equation system (equations ( 3 . 16 ) and ( 3 . 17)) then is 
transformed into the 3 equations 


c + u |a + „ la') 

\dt cbc <3y J 

p f|l + u ^ + v 
\dt cbc dy/ 


pX _ ^ + &\ 

^ \dx 2 dy 2 / 

p y_^ +ti ^ + &\ 

* W * 2 / 


(4.2) 


du dv 
dx dy 


0 


with the three unknown factors u, v, p (X and Y are the components 
of the mass force K per unit volume). 

After various minor transformations the equation system may be 
written as a single equation. To this end one introduces the rotational 
vector rot w which for the plane case has only one component not 
equalling zero: 


- rot z w = co z 



(4.3) 
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Furthermore, the mass force in equation (4.2) is put equal to zero. 
This is permissible in all cases where the fluid is homogeneous and no 
free surfaces are present. In order to introduce o> z into equation (4.2), 
the first equation of (4.2) is differentiated with respect to y, and 
the second with respect to x; then the first is subtracted from the 
second and one obtains: 



With this transformation the pressure terms have been eliminated. 
Equation (4.5) may now, with p/p = V, be written: 

— = V&d (vorticity transport equation) (4.6) 


with cd = o z being denoted as the vorticity. 

This equation signifies: The convective (substantial) variation 

of the vortex strength equals the dissipation of vorticity by 
friction. 

Equation (4.6) forms with the equation of continuity a system of 
two equations with two unknowns, namely u and v, the derivatives of 
which define cd. 

By introducing a flow function t(x,y) one may finally introduce 
a single equation with the unknown i|r. The flow function represents 
the integral of the equation of continuity. One sets: 


That is, therefore, the 
equation of continuity 
is identically 
satisfied by 'if. 
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Moreover, 


a> 




(4.8) 


That is: The Laplacian of the flow function is exactly minus two times 

as large as the vorticity (angular velocity). With this result 
equation (4.5) "becomes, after division "by P: 


(& + u £ + v = VAa ' 1, 


( 4 . 9 ) 


or expressed only in i|r with equation (4.7): 


c&fr + dt c&t _ cW c&>!> 
3t dy dx 8x 8y 


VAA-ty 


(4.10) 


This one equation with the unknown is the vorticity tr an sport 
equation, "but written in terms of ''If. 

The inertia terms are again on the left, the viscosity terms on the 
right side. Equation (4.10) is a differential equation of the fourth 
order for the flow function. Again, its general solution is extremely 
difficult because of the non-linearity. For very slow (creeping) motions 
the friction terms very strongly predominate. Then one may set: 


AAi = 0 


(4.11) 


This simplification is permissible only because the differential equation 
remains of the fourth order, so that no boundary condition is lost. 
However, being linear, this equation is at least solvable. It appears 
alsb in the theory of elasticity where it is designated as the bipotential 
equation. There exists a solution of equation (4.11) by Stokes for 
moving droplets^ which was extended by Cunningham to very small drop 
diameters (comparable to the mean free path of the molecules). 

Herewith we shall conclude the more general considerations and turn 
to the boundary layer problem proper, limiting ourselves to fluids of 
very small viscosity V. 

A few preparatory considerations will lead up to the boundary layer 
problem. One might conceive the notion of simply eliminating all the 
friction terms of the Navi er— Stokes’ differential equation in the case of 
smal 1 viscosity. However, this would be fundamentally wrong as will be 
proved below. 
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An equation which is completely analogous to equation (4.5) occurs 
in the theory of heat transfer: 


c 


P 



+ u 




= X 




(4.12) 


where the velocity components are retained whereas the rotation cn 
replaces the temperature d, the density p the specific heat Cp 

per unit volume, and the viscosity p the thermal conductivity X. On 
the left of equation (4.12) stands the temperature change due to 
convection, on the right the change due to heat transfer. 

The temperature distribution around a heated body immersed in a 
flow with the free stream velocity u Q (for instance fig. 8) is 

determined by the differential equation (4.12). One perceives intuitively 
that for small u Q the temperature increase starting from the body 
extends toward the front and all sides far into the flow (solid contour) 
whereas for large u Q this influence is mainly limited to a thin layer 
and a narrow wake (dashed contour). 

The analogy of equations (4.12) and (4.5) indicates that the 
friction— rotation distribution in question must be similar: For small 

free stream velocity the rotation is noticeable at large distance from 
the body, whereas for large u Q the rotation is limited to the immediate 
neighborhood of the body. 

Thus for rapid motions, that is, large Reynolds numbers (compare 
next section), one expects the following solution of Navi er— Stokes’ 
differential equations: 

1. In the region outside of a thin boundary layer m = 0, that is, 
potential flow 

2. Inside this thin boundary layer cd / 0, thus no potential flow. 

Therefore, one must not set cu = 0 in this boundary layer, even 
for small viscosity. 

It is true that the potential flow is also a solution of Navier— 
Stokes’ differential equations, but it does not satisfy the boundary 
layer condition v^ = 0. 

Proof: The potential flow may be derived from potential $(x,y,z) 

as: 


w = grad $, with A4 = 





(4.13) 
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However, if Al> = 0, then also grad = A gradO = 0, that is. Aw = 0 
for potential flow. According to equation (3.18) this fact signifies 
that in the Navier— Stokes differential equations the friction terms 
vanish identically, and hence that the potential flow actually satisfies 
the Navier-Stokes differential equations. However, it satisfies only 
the one boundary condition v n = 0. 

Thus, for the limiting case of small viscosity, one obtains useful 

solutions for the limiting process v ^ 0 not by cancelling the 

friction terms in the differential equation, since this reduces its 
order (the differential equation of the fourth order for the flow 
function would turn into an equation of the second order; the Navier- 
Stokes differential equations would change from the second to the first 
order), so that one can satisfy only correspondingly fewer boundary 
conditions. 

Thus the limiting process V :> 0 must not be performed in the 

differential equation itself, but only in its solution. 


This can be clearly demonstrated on an example (referred to for 
comparison by Prandtl) of the solution of an ordinary differential 
equation. Consider the damped oscillation of a mass point. The 
differential equation 


m + k + cx = 0 
dt 2 


( 4 . 14 ) 


applies in which m represents the oscillating mass, k the damping 
constant and c the spring constant, (x = elongation, t = time). 

Let for instance the two initial conditions be: 


t = 0; x = 0; dx/dt = 1 


In analogy to the case in question one considers here the limiting 
case of a very small mass m, since then the term of the highest order 
tends toward zero. If one would simply put m = 0, one would treat 
nothing but the differential equation 


k2f+cx = 0 ( 4 . 15 ) 


which by assuming the solution to be of the form x = A e^ is 
transformed into k A, + c = 0, whence A, = -c/k. That is, the 
solution reads: 
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x 



(4.16) 


However, the two initial conditions x = 0 and dx/dt = 1 at the 
time t = 0 cannot he satisfied with this solution. But if one treats 
the complete differential equation (4.14) in the same manner there results 

p 

mX + k X + c = 0 


and hence: 



- k + k \ 1 


2m 


4cm 

k 2 


or the square root might he developed into a series and (since now the 

limiting process m >0 is to he performed) broken off after the 

second term: 




c. 

k 


Thus X-^ corresponds to the previous solution of the first order 
differential equation, where, however, Xg had been lost. For very 
small m, X^ ~ -k/m; therewith the general solution becomes, by 
combination of the particular solutions. 


x 


*1 


X-| t 
e x + 


*2 


e 


= A i 




A 2 e 


m 


(4.17) 


Since for t = 0, x is also supposed to equal zero, there follows: 
A^ = - A^, thus: 


x 



(4.18) 


This equation is plotted schematically in figure 9« The first term of 
equation (4.18), which alone cannot satisfy the boundary conditions. 
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starts from the value A^ at the time t = 0 and decreases expo- 
nentially. The second term is important only for small t— values and 
plays no role for large t. It is very rapidly variable and assures 
that the total solution (solid line) satisfies the boundary conditions. 
The slowly variable solution (in corresponds to the potential flow, 

the second, rapidly variable particular solution (in X^) indicates, as 

it were, the narrow region of the boundary layer; the smaller m, the 
narrower this region. 

Herewith we shall conclude the general remarks and turn to the law 
of similarity. 


CHAPTER V. REYNOLDS’ LAW of SIMILARITY 


So far no general methods for the solution of the Navier-Stokes 
differential equations are known. Solutions that are valid for all 
values of the viscosity are so far known only for a very few special 
cases (for instance, Poiseuille's pipe flow). Meanwhile the problem 
of flow in a viscous fluid has been tackled by starting from the limits, 
that is, one haB treated on the one hand flows of very great viscosity, 
on the other hand flows of very small viscosity, since one obtains in 
this manner certain mathematical simplifications. However, starting 
from these limiting cases one cannot possibly interpolate for flows of 
average viscosity. 

The theoretical treatment of the limiting cases of very great and 
very small viscosity is mathematically still very difficult. Thus 
research on viscous fluids was undertaken largely from the experimental 
side. The Navier-Stokes differential equations offer very useful indi- 
cations, which permit a considerable reduction of the volume of experi- 
mental investigation. The rules in question are the so-called laws of 
similarity. 

The problem is: Under what conditions are the forms of flows of 

any liquids or gases around geometrically similarly shaped bodies them- 
selves geometrically similar? Such flows are called mechanically similar. 

Consider for instance the flows of two different fluids of different 
velocities around two spheres of different size (fig. 10). Under what 
conditions are the flows geometrically similar to each other? Obviously 
this is the case when at points of similar position in the two flow 
patterns the forces acting on volume elements at these points have the 
same ratio. Depending on what kinds of forces are in effect, various 
laws of similarity will result from this requirement. 

Most important for this investigation is the case where all forces 
except the inertia and friction forces are negligible. Furthermore, no 
free surfaces are to be present, bo that the effect of gravity is 
compensated by the hydrostatic pressure. In this case the flow around 
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the two spheres is geometrically similar when the inertia and friction 
forces have the same ratio at every point. 


The expressions for the inertia and viscosity forces acting on the 
volume element will now he derived: there is as friction force per 


unit volume 


4T 

dy 


P u 


du 

dx 


The ratio 



whereas the inertia force per unit volume is 


inertia force 


p u 


friction force 



(5.1) 


must, therefore, he the same at all points of the flow. One now inquires 
as to the variation of these forces with variation in the quantities 
characteristic of the phenomenon: free stream velocity V, diameter d, 

density p, and viscosity p. For variation of Y and d the indi- 
vidual quantities in equation (5.1) at similarly located points vary 
as follows: 


u 




du T 5 2 u 



Therewith equation (5»l) "becomes: 


inertia force 
friction force 



(5.2) 


The law of mechanical similarity is therefore: The flows around 

geometrically similar bodies similarly located and alined with respect 
to the flow have, for equal p Y d/p, geometrically similar stream 
lines as well. If the flows in question are, for instance, two flows 
of the same fluid of equal temperature and density (p and p equal) around 
two spheres, one of which has a diameter twice that of the other, the flows 
are geometrically similar provided that the free stream velocity for the 
larger sphere has half the magnitude of that for the smaller sphere. 
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The quantity p V d/p is, as a quotient of two forces, a dimension- 
less number. This fact is immediately recognized by substituting for 
the quantities their dimensions: 


P 


kg sec 


2 



kg sec 
m 


pVd _ kg sec 2 m m 2 
p m — 4 sec kg sec 


= 1 


This law of similarity was discovered by Osborne Reynolds in his 
studies of fluid flows in a pipe. The dimensionless quantity is called 
after him: 


pV d/p = V d/v = Re = Reynolds* number 


The introduction of this dimensionless quantity helped greatly in 
advancing the development of modern hydrodynamics. 


Connection between Similarity and 
Dimensional Considerations 

As is known, all physical laws can be expressed in a form free of 
the units of measure. Thus the similarity consideration may be replaced 
by a dimensional analysis. The following quantities appearing 
in the Navier— Stokes differential equations are essential for the stream 
line pattern: V, d, p, p. The question is whether there is a combination 

V 1 d P p? p 8 

which is a Reynolds number and therefore has the dimension 1. This 
amounts to determining a, 0, 7, 6 in such a manner that 

v“ d p p* 4° - K° L° T° = 1 (5.3) 

with K, L, T representing the symbols for force, length, and time, 
respectively. Without limiting the generality a may be set equal to 
unity (a = 1) since any power of a dimensionless quantity is still a 
pure number. With a = 1 there results from equation (5«3) 
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= K° L° T° 


(5-M 


By equating the exponents of L, T, K on the left and right sides one 
obtains the three equations: 




K: 

7 + 

8=0 



L: 1+3 

-ky - 

o 

II 

60 

OJ 

> 

(5-5) 

T: 

2 y + 

8 = 1 






The solution gives: 

3 = lj 7=1; 6 - -1 (5.6) 


Accordingly the only possible dimensionless combination of V, d, p, p, 
is the quotient 


^£ = Re (5-7) 

I i 


This dimensional analysis lacks the pictorial quality of the similarity 
consideration; however, it offers the advantage of applicability even 
when knowledge of the exact equation of motion is still missing, if 
there is only known what physical quantities determine the phenomenon. 


CHAPTER VI. EXACT SOLUTIONS OF THE NAVIER— STOKES EQUATIONS 


In general, the problem of finding exact solutions of the Navier- 
Stokes differential equations encounters insurmountable difficulties, 
particularly because of the non-linearity of these equations which 
prohibits application of the principle of superposition. Nevertheless 
one can give exact solutions for a few special cases, mostly, when the 
second power terms vanish automatically. A few of these exact solutions 
will be treated here. 

One investigates first layer flows in general, , that is, flows 
where only one velocity component exists which, moreover, is not 
dependent on the analagous position coordinate, whereas the two other 
velocity components vanish identically; thus for instance: 
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u = u(y, z, t) 
v = w s 0 


( 6 . 1 ) 


The Navier-Stokes differential equations ( 3 . l6 ) are thereby transformed 
into : 

\ 

du ^ „ 3 d , f d 2 u d 2 u 

* 41 W 2 dz 2 


0 = p Y - ^ 
0 7 


0 = pz-| 


> 


( 6 . 2 ) 


while the continuity equation is identically satisfied. 


a. Pipe Flow, Steady and Starting 


1. Steady pipe flow. 


For the case of the pipe lying horizontally the mass forceB are 
everywhere constant, and the equation system (equation (6.2)) yields: 





( 6 . 3 ) 


with the solution 


p = ax + by + cz 


Thus the pressure is a linear function of the position. 

In the case of the pipe standing vertically (fig. 11) the mass 
force is constant in the y— and z— direction and increases in the 
x— direction corresponding to the hydrostatic pressure of the field of 
gravity. Moreover, if one puts Y = 0 and Z = 0; 


dp/dy = dp/dz =0: v = w = 0 
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Thus there remains from equation (6.2) 


u n Sp 

— - = 0 or — = constant 

ax 2 ^ 


The system (equation (6.2)) uhen becomes, under the further assumption 
of steady flow: 


dy 2 



1 

V- dx 


constant 


(6.4) 


One tries a solution of the form 


u = u (a^+ by2 + cz^) 
max 


(6.5) 


with Ujnax representing the velocity at the pipe center. With the 
condition u = 0 at the pipe wall y^ + z^ = r^ there results: 


a = 1; 



and therefore 


with 


u(y, z) 



r 2 dp 

U »ax ~ dx 



( 6 . 6 ) 


This solution is identical with equations (2.3) and (2.4) if one bears 
in mind that in those equations y represented the radial distance from 
the pipe center. Thus Poiseuille's pipe flow was found as an exact 
solution of the Navi er— Stoke s differential equations. 


2. Starting flow t h rough a pipe . 

By the expression, "Starting flow through a pipe," the following 
problem is meant: Let the fluid in a circular pipe of infinite length 
be at rest until the time t = 0. At the time t = 0 let a pressure 
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difference invariable with time, - p^, suddenly be established 

across a piece of pipe of length 2, so that for the entire pipe the 
following equation applies: 


t > 0: 


6p 

dx 



= constant 


( 6 . 7 ) 


Here 



and the inertia and friction terms balance each other. 


The solution of this problem, which will not be further discussed here, 
was given by Szymanski (reference 12). The phenomenon is not dependent 
on the longitudinal coordinate x. The velocity profiles at various times 
can be seen from fig. 12. It is characteristic that, first, at the 
pipe center, the velocity remains locally constant and the friction is 
noticeable only in a thin layer near the surface. Only later does the 
friction effect reach the pipe center. Poiseuille*s parabolic profile 
of the steady pipe flow is attained asymptotically for t-> » . 


One must clearly distinguish between the non— steady starting flow 
through a pipe discussed here and the Bteady pipe inlet flow. This 
latter is the flow at the inlet of a pipe. The rectangular velocity 
profile present in the entrance cross section is, with increasing 
distance x from the inlet, gradually transformed under the effect of 

friction into Poiseuille’s parabolic profile. Since here ^ 4 0, 

ox 

this is not a layer flow. This inlet flow was, for the plane problem, 
exactly calculated from the differential equations by H. Schlichting 
(reference 14) and for the rotationally-symmetrical problem, according 
to an approximate method, by L. Schiller (reference 32). 


Third lecture (Dec. 15, 1941) 
b. Plane Surface; a Surface Suddenly Set in Motion 
and an Oscillating Surface 
1. Oscillating Surface. 


Let a plane surface of infinite extent perform in its plane recti- 
linear oscillations in the x— direction (fig. 13). The y-direction is 
assumed normal to the surface. Let the oscillation take place with the 
velocity u Q = A. cos n t, with A denoting the amplitude, n the 

frequency of the oscillation. The fluid near the surface is carried along 
by the friction. Since the surface is of infinite extent, the state of 
flow is independent of x and z. The flow in question is therefore a 
non— steady plane layer flow for which 
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u = u(y,t); 


v = w = 0; — = 0 

8x 


If, finally, one puts the mass forces equal to zero, the Navier-Stokes 
differential equations (equation (3.16)) are reduced to 


da = V (6.8) 

8t 2 

dy 


with the boundary condition y = 0, u = u Q cos n t. This equation has 

the same structure as the differential equation for thermal expansion in 
a rod (linear equation of heat conduction). It becomes identical with 
it if v is replaced by the thermal conductivity a and u by the 
temperature (compare equation (4.12)). 

For the solution of equation (6.8) one uses 

u = Ae - ^ cos (nt — ky) (6.9) 


where k is a constant to be determined. 
Then 


di = — n Ae~^ sin (nt — ky) 
ot 


di = - k Ae ^ cos (nt — ky) + k Ae~^ sin (nt — ky) 

oy 

= — 2k 2 Ae""^ sin (nt — ky) 

dy 2 


Insertion in equation (6.8) gives nA = 2k 2 AV, and thence 


k = \| Ja- 
2 v 


(6.10) 


The velocity distribution u(y, t) is therefo re an oscillation with 


amplitude decreasing toward the outside Ae * 2 * j the layer at the 
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distance y from the surface has a phase lag of 



with respect to 


the motion of the surface. The wave length X in the y— direction 
follows from kX = 2« as 


X 



( 6 . 11 ) 


Thus the co-oscillating layer is thinner, the greater the frequency n 
and the smaller the viscosity V. The result X ~ is to he noted. 

The velocity profiles for various times are given in figure 13 . 

2. The Surface Suddenly Set in Motion, 

Equation (6.8) yields another exact solution of the Navier— Stokes 
differential equations, namely the flow in the neighborhood of a plane 
surface which suddenly starts moving in its own plane with the constant 
velocity u Q . In this case the boundary conditions are 


t — 0, for all y: u = 0 


t > 0, y = 0: u = u 


( 6 . 12 ) 


An appropriate variable for the solution is the dimensionless quantity 


T) = 



( 6 . 13 ) 


As the solution of equation (6.8) with the boundary conditions 
equation (6.12) one obtains 


u(y, t) 



(6.14) 


The correctness of the solution is readily confirmed by substitution. 
The variation of the flow with time is indicated in figure 14. 

The probability integral appearing in equation (6.14) -§r 

V* 

has, for 11 = 7] =1.9, the value 0.01. Therefore, for tj = q , 

o 

u = 0.01 u Q . The thickness of the layer carried along by the friction 
is, therefore. 



dq 
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5 " 2 \ V'vt 


(6.15) 


Thus here again the friction layer thickness is 


6 ~ \j~V 


At this point we conclude the discussion of layer flows and turn to a 
few other exact solutions of Navier-Stokes differential equations. 


c. Plane Stagnation-Point Flow 

The plane flow in the neighborhood of a stagnation point on a smooth 
wall is considered. With the coordinate system according to figure 15 
the corresponding potential flow has the potential 

• = | (x 2 - y 2 ) (6.16) 


and the stream function 


t = a x y 


The velocity components are: 


U = ax; V = -ay (a = constant) 


(6.17) 


This is a potential flow which, coming from the direction of the y— axis, 
encounters the solid wall y = 0, divides, and flows off parallel to 
the x-axis. Whereas the potential flow glides along the wall, the viscous 
flow must adhere to it. If one designates the velocity components of 
the viscous flow by u(x, y) and v(x, y), the boundary conditions for 
them are: 


y = 0; u = 0; v = 0 
y = °°; u = U; v = V 


(6.18) 


For the stream function of the viscous fluid one uses the equation; 
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t (x, y) = x f (y) 


( 6 . 19 ) 


Thence one obtains: 


u = x ^ = x f*(y) 

dy 


v = - f(y) 


( 6 . 20 ) 


The boundary conditions (Equation ( 6 . 18 )) require 


y = 0: r = f* = 0 
y = co : f = ay; f ’ = a 


( 6 . 21 ) 


At large distances from the surface, that is, in the potential flow, the 
pressure is calculated from Bernoulli’s equation 

p = P o - | W 2 = p o - | a 2 (x 2 + y 2 ) (6.22) 

with W = V 2 signifying the magnitude of the velocity of the 

potential flow, given by equation ( 6 . 17 ). Let p Q be the total pressure 
of the potential flow. 

For the viscous flow one formulates the analogous equation 


P = Po - f a 2 [x 2 + F(y)] 


( 6 . 23 ) 


With these equations one turns to the Navier— Stokes equations which read, 
for vanishing mass forces: 


u 



+ 


du 
v — 


dy 



dx 



dy 


1 dp 
P dx 



+ 



i§E + v 

p dy 




! 

J 


( 6 . 24 ) 
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The continuity equation has already been integrated by introduction of 
the stream, function 4'. By substitution of equations (6.20) and (6.23) 
into equation (6.2k) one obtains 

x(f’2 — f f") = a 2 x + Vxf*" 

2 (6.25) 

f f» = — F*(y) - V f" 


These are two differential equations for the two unknown functions 
f(y) and F(y) which determine the velocity distribution and the 
pressure distribution, respectively. The component of the velocity u 
parallel to the surface, that is, the function f*(y), is of particular 
interest. Since F (y) does not appear in the first equation, one 
solves first the first equation and then, after substitution of this 
solution, the second. Thus the differential equation to be solved reads 
at first 


- f f " = a £ 


+ v 


ft** 


(6.26) 


with the boundary conditions according to equation (6.21). 

This non-linear differential equation cannot be solved in closed 
form. If one introduces instead of y the variable 


I = a y 


( 6 . 26a ) 


and in addition the similarity transformation 


f (y) = Atp(5) 


( 6 . 26b ) 


Then the inhomogeneous term in equation (6.26) becomes equal to 1, and 
the solutions therewith become independent of the specific data for the 
flow. Thereby equation (6.26) becomes 


a2 A 2 


/ d® 


- <P 


d 2 qp 

dt 2 


. a 2 + VAa} & 

dS 3 


( 6 . 2 ?) 


pop o p 

If one now equates a A = a and vAa J = a , that is. 
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(6.27a) 


the differential equation for <p(5) with I 



reads : 


<P 


» » 


cp <p" = cp' 2 + 1 = 


(6.28) 


with the "boundary conditions 


5 = 0: cp = <p* = 0 

| = oo; (p* = 1 


The solution found "by series development can "be found in the thesis of 
Hiemenz (reference 10) , compare table 1* . The velocity component 
parallel to the Burface is 


5 = ^ f*(y) = <P'( I) (6.29) 

U a 


It is indicated in figure l6. The curve cp* ( |) increases linearly at 
1 = 0 and approaches one asymptotically. For about | = 2.6, cp’ ~ 0.99; 
thus within about one percent of the final value. If one again designates 
the ‘corresponding distance from the surface y = 6 as the boundary layer 
thickness (friction layer thickness), then 

6 - | ‘o - 2 - 6 \Jz (6.30) 


Thus in this flow, as in the former ones. 


5 


It is also remarkable that the dimensionless velocity distribution 
according to equation (6.29) and the boundary layer thickness according 
to equation (6.3O.) are independent of x, thus do not vary along the wall. 


*The tables appear in appendix, chapter XII. 




38 


NACA TM No. 1217 


For later applications the characteristics important for the friction 
layer, displacement thickness 5* and momentum thickness d, are 

introduced here; they are defined "by 


U5 * = f° (U — u) dy (6.31) 

y=o 


U 2 d = f u(U - u) dy (6.32) 

y=o 


The displacement thickness gives the deflection of the stream lines of 
the potential flow from the surface hy the friction layer; the momentum 
thickness is a measure of the momentum loss in the friction layer. By 
insertion of equation (6.29) in (6*31) and (6.32) and calculation of the 
definite integral one finds 


6 * 


■a 



(1 - cp* ) as = 0.6482 ^ 


9* (1 - cp*) dl = 0.2923 


Vs 


(6.33) 


(6.3^) 


and hence 


**= 2.218 

■a 


(6.35) 


The quantity &* is indicated in figure l6. For comparison with a 
later approximate solution one also notes the numerical value of the 

dimensionless quantity ^7- — . One finds from equations (6.17) and (6.33) 

v dx 

— = 0.4202 (6.36) 

V (Jx 


The exact solution of the Navier— Stokes differential equations found 
here gives, theref ore , • f or large Reynolds numbers a friction layer 

thickness decreasing with V; and a transverse pressure gradient 
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decreasing with pa \/v&* Both confirm the boundary layer assumptions* 
to he discussed later. 


d. Convergent and Divergent Channel 

A further class of exact solutions of the Navi er— Stokes differential 
equations exists for the convergent and divergent channel with plane walls 
(fig. 17), as given hy G. Hamel (reference 11). 

Without entering into the details of the rather complicated calcu- 
lations the character of the solutions will he briefly sketched: 

The velocity distrihutions for convergent channels, plotted against 
distance along the surface for various included angles a. and for various 
Re-numbers appear as indicated in figure 18. At the tunnel center the 
velocity is almost constant, and at the surfaces it suddenly declines to zero. 

In the case of divergent tunnels one obtains greatly differing 
forms for the velocity profiles, depending on the included angle and 
the Re-number. Here all velocity profiles have two inflection points. 

For small Re-nunibers and small included angles the velocity is positive 
over the entire cross section (solid curve in fig. 19); for larger angles 
and larger Re-numbers, on the other hand, the velocity profiles have 
reverse flow at the surface (dashed curve in fig. 19)* The reverse flow 
is the initial phase of a vortex formation and therefore of the separation 
of the flow from the surface. Generally, the separation does not occur 
symmetrically on both surfaces; the flow separates from one side and 
adheres to the other surface (fig. 20). 

These examples also confirm the theory that exact solutions have the same 
character as approximate solutions of boundary layer theory; in particular, 
they confirm that for the convergent channel a very thin layer with con- 
siderable friction effect is present near the surface (here also the 
calculation shows that the layer thickness \fv) and that for the 
divergent channel reverse flow and separation occur. 

We here conclude the chapter on the exact solutions of the Navier— 

Stokes differential equations and turn to the approximate solutions. 

By exact solutions have been understood those where in the Navier-Stokes 
differential equations all terms are taken into consideration that, in 
the various cases, are not identically zero. By approximate solutions 
of the Navier-Stokes differential equations will be understood, in 
contrast, solutions where terms of small magnitude are neglected in the 
differential equations themselves. However, by no means are all the 
friction terms to be neglected simultaneously, since this would represent 
the case of potential flow. 


*The rotationally— symmetrical stagnation— point flow has been 
calculated by Homann (reference 17). Instead of equation (6.28) one 

obtains the differential equation q) ,M + 2qxp" — qp * 2 + 1 = 0. 
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CHAPTER VH. VERY SLOW MOTION (STOKES, OSEEN) 


The exact solutions of the Navi er— Stokes differential equations 
discussed in the previous chapter are of a very special kind. Most of 
them dealt with flows along a plane surface, where the stream lines 
are rectilinear. Most flows existing in practice, as for instance flows 
around arbitrary bodies, cannot be calculated exactly from the Navier— 
Stokes differential equations, but must be treated by approximate methods. 
Two kinds of such approximations are possible: 

1. For predominant viscosity, completely neglecting the inertia terms 
suggests itself (very small Re-number; Re < 1). 


2. For very small viscosity and therefore predominant inertia one 
takes the viscosity into consideration only in a very thin 
layer in the neighborhood of the solid wall; for the rest, 
the flow is regarded as frictionless. Here the Re-number is 
very large (Prandtl’s boundary layer theory). 

The first limiting case with very small Re-Cumber will be discussed 
in this chapter. A small Re-number indicates small velocities, small 
body dimensions, and large viscosity. Since the inertia terms depend 
on the square of the velocity whereas the friction terms are linear, 
all inertia terms in the Navier-Stokes differential equations are, for 
very small Re— numbers, negligible. It is to be expected that an 
approximation will thereby be obtained for very slow (creeping) motion, 
as for instance the falling of a minute fog particle*) or the slow 
motion of a body in a very viscid oil. 


Neglecting all inertia terms one obtains from the Navier-Stokes 
differential equations ( 3 • l6 ) the following: 


dp 

dx 


MAul 


! = mat 


dp 

bz 


MAwJ 


(7.1) 


du dv dw 

dx dy dz 


(7.2) 


*For a sphere falling in air (v = lh X 10 -6 m 2 /sec) for instance: 
Re = Y d/v = 1, for d = 1 mm; V = l.hO cm/sec. 
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The same boundary conditions apply to this system of equations as apply 
to the complete Navier— Stokes differential equations, namely vanishing 
of the normal component v n =■ 0 and the tangential component v^ = 0 

at the bounding surfaces. 

The neglect of all inertia terms in Navier— Stokes differential 
equations does not represent as serious ah inaccuracy as the neglect of 
all friction terms when transforming the Navier-Stokes differential 
equations into Euler’s differential equations of the frictionless flow. 
That is, by neglecting the inertia terms, the order of the differential 
equations is not lowered so that in the simplified differential equations 
the same, boundary conditions as in the Navier-Stokes complete differential 
equations can still be satisfied. 

Furthermore one obtains from the equations (7*1), taking into 
account the continuity, by differentiating the first with respect to x, 
the second with respect to y, the third with respect to z, the following 
equation for the pressure p 

^E + 2!e + ^£ = A P = 0 (7.3) 

<* 2 a/ bz 2 


that is, for creeping motions the pressure function p(x, y, z) is a 
potential function. 

The details of the calculation will not be discussed more thoroughly, 
particularly since the creeping motion is technically not very important. 
However, at least Stokes* famous solution for the sphere will be discussed 
briefly (fig. 21). The drag of a sphere for creeping motion consists of 
the contributions of the pressure drag (form drag) and the surface fric- 
tion drag. The latter is obtained by integration of the wall shearing 
stress over the entire sphere surface. Stokes performed the integration 
of the equation systems (7.1) and (7*2) for a sphere in a uniform 
flow of velocity TJ 0 . There results, according to Stokes, for the entire 

drag of the sphere of radius R: 


W = W_ + W_ = 6 jqjR U (7.^) 

Dr R o 


The drag is* therefore* proportional to the first power of the velocity. 

If one introduces for the sphere a drag coefficient c w which* in the 
customary manner* is referred to the frontal area and the dynamic pressure 
of the free stream velocity 



! 
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W = 


"w 


it K 


£ n 4 

2 u o 


(7.5) 


there results for the drag coefficient according to Stokes’ formula: 

<V--i (7-6) 

T?p V 


One can state immediately that the stream line pattern of this 
creeping motion must he the same ahead of and behind the sphere since 
for reversal of the initial flow (sign reversal of the velocity 
components) the equation system (equation ( 7 . 1 )) goes over into itself. 
The stream— line pattern for the viscous sphere flow, as it presents 
itself to an observer who is at rest relative to the flow at infinity, 
is shown in figure 22. The fluid particles are pushed aside by the 
sphere in front and come together again behind it. 

As shown by a comparison of Stokes' drag formula equation (7-6) 
with test results (reference 33 )j this formula is valid only for the 
region Re < 1. 


Correction by Oseen 

In Oseen 1 s later improvement of Stokes' solution for the sphere 
.the inertia terms in the differential equations are partly taken into 
consideration. Oseen formulates the velocity components u, v, w; 


u = U 0 + u' ; v = v'; w = w' 


(7.7) 


where u' , v' , w* may be considered as disturbance velocities which in 
general are small compared with the free stream velocity U Q . This 

assumption is not actually correct for the immediate proximity of the 
sphere surface. With the formulation (equation (7»7)) the inertia terms 
in equation ( 3 . 16 ) are divided into two groups, for instance: 



dv 

0 ST' 


, t du 1 t dv 1 

and , u f , . . . 

Sx &c 


The second group of second order, as compared with the first group, is 
neglected. Therewith one then obtains from the Navier— Stokes equations 
of motion ( 3 . 16 ) the following equations of motion, which are taken as 
a basis by Oseon. 
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p U 0 -^ + ^ = pAu’ 
OX ox 


TT ^ yl . At l 

P U 0 - 7 — + ^ ■ uAV» > 


boa by 


.. dw' dp . , 

p IT + — = pAw' 

0 dr dz 


In addition, one usee the continuity equation: 


du' dv* + dw* _ q 

dx by dz 


(7.8) 


(7.9) 


and the same "boundary conditions as in the Navier-Stokes differential 
equations. One calls the contributions of the convective terms in these 
equations that were taken into consideration, for instance U 0 the 

dx 


semi— quadratic terms. These differential equations of Oseen and Stokes' 
differential equations are both linear. The stream line pattern, as it 
results for this sphere flow according to Oseen, is given in figure 23 • 
Here again the observer is at rest relative to the fluid at large distance 
from the sphere. Thus the sphere is dragged past the observer with the 
velocity Uq. The stream line pattern ahead of and behind the sphere are 


now not the same, as was the case in Stokes' solution. Ahead of the 
sphere exists almost the same displacement flow as in Stokes' pattern; 
behind the sphere, however, the stream lines are closer together, that is 
the velocity is greater here than in Stokes' case. A wake is present 
behind the sphere similar to that from test results for large Reynolds 
numbers. 


For the sphere drag calculated by Stokes there results with the 
drag coefficient c w introduced in equation (7.5) "the formula: 


c„ = ft (l * 

Re ^ 


J-Re^jRe ' 

16 ) V 


( 7 . 10 ) 


The test results (reference 33 ) show that Oseen' s formula is fairly 
accurate up to about Re = 5« 

With these brief re mar ks we conclude the limiting case of small 
Reynolds numbers and turn to the case which is of foremost interest in 
practice: the case of very large Reynolds number. 
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CHAPTER VIII. PRANDTL’ S BOUNDARY LAYER EQUATIONS 


The other extreme case of very email viscosity or of very large 
Reynolds number will now he treated. In this case the inertia effects 
are predominant within the main body of the fluid whereas the viscosity 
effects there are almost negligible. 

A significant advance in the treatment of motion of fluids for 
large Reynolds numbers, that is, in general, of fluids of very small 
viscosity, was attained by L. Prandtl in 190k (reference 7). Prandtl 
demonstrated in what way viscosity is essential for large Reynolds 
numbers and how one can simplify the Navier— Stokes differential equations 
in order to obtain at least approximate solutions. 

Let us consider the motion of fluid of very small viscosity, for 
instance of air or water surrounding a cylindrical streamline body 
(fig. 2k). Up to very near the surface the velocities are of the order 
of magnitude of the free stream velocity U Q . The stream line pattern 

as well as the velocity distribution agree to a large extent with those 
of the frictionless fluid (potential flow). More thorough investigations 
show, however, that the fluid by no means glideB along the surface (as 
in potential flow) but adheres to it. The transition from zero velocity 
at the surface to the fully developed velocity as it exists at some 
distance from the body, is effected in a very thin layer. Thus one must 
distinguish between two regions which, it is true, cannot be rigorously 
separated: 

1. A thin layer in the immediate proximity of the body where the 

velocity gradient normal to the surface “ is very large 

&n 

( boundary layer ). Here the viscosity p, though very small, 
plays an essential role inasmuch as the frictional shearing 

stress t - (i — can assume considerable values, 
on 

2. In the remaining region outside of this layer velocity gradients 

of such magnitude do not occur, so that there the effect of 
viscosity becomes insignificant. Here frictionless potential 
flow prevails. 

In general 'one may say that the boundary layer is thinner, the 
smaller the viscosity or, more generally, the larger the Re-number. It 
was shown before on the basis of exact solutions of the Navier— Stokes 
differential equations that the boundary layer thickness is 

5 ~ \J v 

The approximations to the Navier— Stokes differential equations to be made 
below are more valid the thinner the boundary layer. Thus the solutions 
of the boundary layer equations have an asymptotic character for infinitely 
increasing Reynolds numbers. 
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Let us now make the simplifications of the Navier-Stokes differential 
equations for the boundary layer . To this end the order of magnitude of 
the separate terms of the Navier-Stokes differential equations must he 
estimated. One considers the flow around a cylindrical hody according 
to figure 24. One imagines the Navier-Stokes differential equations 
written non-dimensionally, hy referring all velocities to the free stream 
velocity U 0 and the lengths to a hody length 2. The pressure will he 
made dimensionless with p U 0 2 , the time with l fU Q . Furthermore 

Re = —2— represents the Reynolds number. Accordingly, the Navier-Stokes 
V 

differential equations become — omitting the mass forces according to 
equation ( k.2 ), hy writing the same letters for the dimensionless quantities 
as for the dimensional ones — 


Su Su Su 

— + u — + v — 

St Sx Sy 

1 11 5 ^ 


Sv Sv Sv 

— + u — + v — 

St Sx Sy 

6 15 5 1 

Su | Sv 
Sx Sy 

1 1 

The estimation gives: Longitidinal velocity u is of the order of 

magnitude 1. Dimensionless houndary layer thickness S/2«l. Therefrom 
follows : 

Su _ 1. S 2 u ^ _1_ 
ay 2 S 2 

whereas the derivatives with respect to x are of normal order of 
magnitude, thus 

Q 

Su . S' u 

— ~ 1 j — ~ x 

31 a ^ 2 
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From the continuity equation follows therefore: 




& ~ 8 
8y 


The transverse velocity v in the boundary layer is therefore, to the 
first order, small in comparison with the longitudinal velocity. Further 
there follows: 


— 

chc 


~ 6; 


d 2 V 

— 2 ~ 6j 

afv 

2 


dy 


Thus there result for the single terms of Navi er— Stokes differential 
equations the orders of magnitude noted in equation (8.1): On the right 

side in the first equation 




so that it can be neglected 

One now obtains within the boundary layer friction terms which are 
of the same order of magnitude as the inertia terms, if l/R is of the 
order of magnitude 5^, or, if the dimensional quantities are again 
written, 



( 8 . 2 ) 


In the second equation of motion all terms then are of the order of 

rha 

magnitude 5, including the transverse pressure gradient — . In the 

8y 

boundary layer, as long as it is thin, the dependence of the pressure 
on y may therefore be neglected. Thus approximately the same pressure 
prevails within the boundary layer as at its edge, that is, the pressure 
of the potential flow. The pressure w ithin the boundary layer is therefore, 
as it were, impressed by the potential flow. 


The second equation of motion is therewith exhausted and does not 
have to be considered further. Using again the dimensional quantities 
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one now makes the Navier— Stokes differential equations assume the following 
simplified form: 


du 8u du 1 dp ,, d^u 

— + u — + v — + v - 

St Sx Sy p ox ^y2 


Su | Sv 
Sx Sy 


0 


y 


Boundary conditions: y = 0: u = v = 0 

y =oo : u = U 


(8.3) 


The pressure in the Boundary layer, which is dependent only on x, is 
determined from the potential flow U (x, t), assumed to he known 
according to Bernoulli* s equation. 


or 


P + f ^ - Pc 


+ £02 


I = U — + — 

p Sx Sx St 


for stationary flow 'l 


> (8.4) 


for nonstationary flow 


/ 


With the potential flow known, the equation system (8.3) represents 
a system of two equations with the two unknowns u and v. 

Numerical example : In order to help clarify the concepts a numerical 

example is given for the thickness of the boundary layer. The problem is: 
What is the boundary layer thickness, for instance, in the caBe of the 
plate in longitudinal flow at the distance 2 = 100 centimeters from the 
leading edge? Let the velocity be Uq = 30 meters per second and the 

kinematic viscosity for air V = 0.14 x lO - ^ meters square per second; 

then the Reynolds number is R = U 0 ^- = 2.1 x 10 and Jr = 1.4-5 x 10 . 

A numerical factor is still missing in the formula (8.2) for the boundary 
layer thickness. For the plate in longitudinal flow it is, as later 
calculations will show, five, provided one understands by the boundary 
thickness 8 that distance from the surface where the velocity has the 
value 0.99 u 0 . Thus a calculation by the formula 

8 5.0 

1 n/R 
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results for the present case in a "boundary layer thickness & = 3*45 “• 

It should he added that the Reynolds number is already so large that the 
boundary layer at the end of the plate would he turbulent. The transition 
from laminar to turbulent lies further upstream, and at this point the 
boundary layer thickness would then be somewhat smaller than the value 
determined above. 


Fourth lecture on December 22, 1942 


Physical S ummar y and Conclusions 


The physical content of the considerations so far can be summed up 
in the following sentences: 

1. In a very thin layer on the body, the boundary layer, the 

velocity passes from the value zero at the surface to the 
value which the potential flow would have in the neighborhood 
of the surface. 

2. The pressure in the boundary layer is practically independent 

of the coordinate normal to the surface and equals the 
pressure of the potential flow along the surface. 

3. In the boundary layer the only friction force to be taken into 

consideration is the shearing stress t = p. — . 

oy 

4. (without proof) The curvature of the surface may be neglected 

in the boundary layer as long as the radius of curvature is 
large compared with the boundary layer thickness (Boltz, 
Thesis, (reference 9)). 

5. All these considerations are valid only as long as no separation 

of the flow from the surface occurs. 

Without integration of the boundary layer equations one can draw from 
these sentences important physical conclusions as to the flow pattern: 

In particular separation occurs if a transport of boundary layer material 
into the interior of the fluid takes place. If a region with pressure 
increase exists along the body contour, the retarded fluid in the boundary 
layer is in general, because of its small kinetic energy, not able to 
penetrate too far into the region of higher pressure. It then withdraws 
laterally from the region of higher pressure, separating from the body, 
and is deflected into the interior (fig. 25). As the point of separation 
one defines the boundary between forward flow and reverse flow of the 
layer nearest the surface, thus 
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point of separation: 



(8.5) 


Determining when and where separation occutb requires for each case the 
integration of the "boundary layer equations. 


One can readily understand that for the velocity profile u(y) at 
the separation point and for all velocity profiles in the decelerating 

> o) an inflection point** must he present. From equation (8.3) 

namely, for the surface y = 0 there follows immediately the relation 


flow 


( 





dp 

die 


( 8 . 6 ) 


The curvature of the velocity profile at the surface therefore exchanges 
signs with the pressure gradient. Thus for flow with pressure decrease 


[ accelerated flow 


V 

also 


8 2 u 


<•>(§). 


<0 is valid and therefore 


dy c 


surface 

<0 in the entire boundary layer (fig. 26), 


For the region 
> 0. 

surface 

However, since in any case at larger distances from the surface 


Of the"preseure Increase (decelerate flow, St > o) (0) 


2-Ji < 0, there must exist* for decelerating flow* within the boundary 

Sy2 a 2 u 

layer a point where — - = 0 (inflection point) (fig. 27). For 

dy 

decelerating potential flow the boundary layer profile has* therefore* 
an inflection point. Since the separation profile with vanishing surface 
tangent must necessarily have an inflection point* it follows that 
separation can occur only when the potential flow is decelerating 


(£ > 0 )- 


* The velocity profile at the point of separation therefore starts 
with a vanishing tangent (fig. 25). Velocity profiles behind the point 
of separation have reverse flow in the neighborhood of the surface (fig. 25). 

** The presence of an inflection point is significant for the 
stability of the velocity profile (transition from laminar to turbulent, 
compare chapter XXI Part II.) 
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If separation is present, the potential flow can no longer envelop 
the body closely everywhere. Thus the pressure distribution sometimes 
deviates considerably from that given by the potential theory. In such 
cases the pressure variation impressed on. the boundary layer can in 
most cases be determined only empirically, because the frictionless outer 
flow itself depends on the complicated phenomena connected with the 
separation. 

Thus the boundary layer theory explains also the fact that in 
addition to the frictional drag a pure pressure drag, called "form drag,” 
appears. 

In regard to later calculations the following explanation shall be 
given: If equation (8.3) is differentiated with respect to y, there 

results for stationary flow 


2 2 2 ? 

o u du c3u o u 8v cfru _ _ 1 op ,, cru 

dxdy dy dx 8y dy p dxdy + ^3 


(8.7) 


Due to ^- = 0 and 
dy 

from equation (8.7) 


to the boundary 
for the surface 


conditions u = v = 0 
y = 0 the relation 


one obtains 



( 8 . 8 ) 


which 1 b valid for all stationary boundary layer profiles (pressure 
increase and pressure decrease). 


Frictional Drag 

As a result of the integration of the boundary layer equation one 
obtains the velocity distribution and the separation point and can there- 
from calculate the particularly interesting surface friction drag in the 
following manner. The friction drag Wg results from the integration of 
the surface shearing stress over the surface of the body. For the plane 
case one obtains, with the symbols according to figure 28, for the friction 
drag 
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The integration has to he extended over "both sides of the surface from 
the stagnation point to the separation point (s = 2*), h signifies the 

height of the cylindrical body. Because cos q> ds = dx and t = 4 

° W / 0 

one obtains for the friction drag 


W R = 2bu 



dx 


( 8 . 10 ) 


tx=o 


This integration alBO haB to be extended along both sides of the body. 
For calculation of the friction drag one needs, therefore, the velocity 
gradient at the surface. The latter can only be obtained by integration 
of the boundary— layer differential equations. If the separation point 
appears ahead of the trailing edge, the formula has to be applied only 
up to the separation point or, sometimes, up to the point of laminar— to— 
turbulent transition which is located further upstream. Behind this 
transition point exists turbulent surface friction drag. In order to 
obtain the total drag, the form drag haB to be added to this friction 
drag; however, the form drag cannot be obtained from the boundary layer 
calculation in a simple manner. 


CHAPTER IX. EXACT SOLUTIONS OF THE BOUNDARY LAYER 
EQUATIONS FOR THE PLANE PROBLEM 
a. The Flat Plate In Longitudinal Flow 


One of the simplest examples for the application of the boundary 
layer equations ( 8 . 3 ) is the flow along a flat plate. Let the plate 
begin at x = 0, extend parallel to the x— axis, and be infinitely 
long (fig. 29). Let the stationary flow of the free stream velocity U D 
be treated. The calculations for it were made by H. Blasius (reference 8) 
in his Gottingen thesis. In this case the velocity of the potential flow 
is constant, thus dp/dx = 0. The boundary layer equations ( 8 . 3 ) become 
therefore 


du . _ 8u _ y 

dx dy 2 

ay 

= 0 

dx dy 


*\ 

> 


✓ 


(9.1) 
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with the boundary conditions 


y=0: u = v = 0 


y = oo : u = U Q 


(9.1a) 


Since the entire system has no characteristic length, the assumption 
suggests itself that the velocity profiles at various distances x are 
affine to each other, so that one may write u/tJ 0 = cp (y/6). Let therein 
5 = 6(x) represent a measure of the boundary layer thickness, increasing 
with the length of run. One arrives at an estimate of the boundary layer 
thickness in the following manner: 

According to the former exact solutions of Navier— Stokes equations 
(chapter VI), for instance for the non— stationary problem of the surface 
suddenly set in motion. 



with t denoting the time since the beginning of the motion. Applied 
to the present stationary problem one may substitute for the time t the 
time required by a fluid particle to travel from the leading edge of the 

plate to the point x. This time t = — and one has, therefore, for 

u o 


the present case 


S~ 



(9.2) 


Thus it is useful to introduce for the distance from the surface y the 
new dimensionless coordinate tj = y/6 or according to equation (9*2) 


T| = y 



(9-3) 


For further calculations one observes that 


= _ i a. = 

dx 2 X* dy 



The continuity equation is again integrated by introduction of the stream, 
function \(r; for the latter one assumes 
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if = fix U 0 f(q) 


(9-4) 


Accordingly f(q) is a dimensionless stream function; for the velocity 
components one obtains from equation (9.4) 


St d^ Sq , v 

u = * * ^ ' u o f (i) 


(9-5) 


7 - - S - - \P~ 5 o f'(,) fa - 




V u, 


2 (Tjf* - f) 


Furthermore one obtains 


(9-6) 


Su _ 1 u 0 -tt. Su tt tl U o ^tt. S^U w o 


U„ 


UU _ -L U «ft # UU TT ^ _p tf , V ^ ^ «,!* 

s7 _u °|fc ' ^2 _ vx f 


(9-7) 


By insertion of these values in equation (9 . 1 ) there results 

U 2 JJ 2 U 

— riff" + -2- (qf *— f ) f" =V-£ f»" 

2x 2x 1 xv 


or, after simplifying, the following differential equation for the stream 
function f(q) 


f f" + 2 f»" = 0 


(9-8) 


Because of equation (9.1a) and of equations (9*5) and (9*6) the boundary 
conditions are 


q = 0: f = f* = 0; q = °°: f ' = 1 


(9-9) 


For the present case therefore, there results from the two partial 
differential equations (9.1), by the similarity transformation (equation 
(9*3) )> 311 ordinary non-linear differential equation of the third order. 
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The three "boundary conditions according to equation ( 9 . 9 ) are sufficient 
for a unique determination of the function f(q) from this equation 

A particular solution of the differential equation ( 9 . 8 ) is the 
solution 


f = T) + constant 


This solution corresponds to the potential flow; we shall revert to it 
presently. 


The general solution of the differential equation (9*8) cannot he 
given in closed form. Therefore one must calculate it either by numerical 
methods or by series developments. Blasius obtained the solution by a 
power series development near T) = 0 and an asymptotic development 
near q = “, which are combined at an appropriate point. Since the 
method of calculation is characteristic of the solution of the boundary 
layer differential equations it will be described in more detail. The 
power series around q = 0 is formulated in the form 

S. 

A A 

f(T)) = A q + A 1 T) + T) 2 + ^ Tl 3 + • . • 


f"(Ti) = A 2 + A 3 T) + ^ T)2 + ^ n 3 + 
f»"(Ti) = A + A^t) + ^ H 2 + “ *1 3 + 

J 4 2 ! 3 ! 

Because of the boundary conditions for q = 0 one 


> 


/ 

has immediately 


( 9 . 10 ) 


A = A = 0 
o 1 


By insertion of the equations (9-10) into the differential equation (9*8) 
one obtains 

2a 3 + n 2A4 + ir ( a 2 2 + 2A 5) + { kA 2 a 3 + 2a 6 ) + . . . = 0 

In order to make the equation ( 9 . 10) represent a solution of the differ- 
ential equation it is required that in the last equation all coefficients 
of the single powers of t \ vanish. First, one has immediately = 0* 

and further 

*This also follows at once from equations (8.6) and (8.8). 
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Aj f A 2 > A 6 - °> Ay - 0 

'T ^ 3 


Thus only the coefficients Ag, A^, Ag. . . . are different from zero. 

The coefficient A 2 remains at first undetermined since the third boundary 
condition for tj = » was not yet satisfied. The remaining coefficients 

A , A_, A . . . can all be expressed in terms of A 0 . One therefore 
5 8 11 d 
sets up, with = <x> a. series for f(ij) which progresses by powers 

of T)^, in the following form: 


OO 



11=0 


(9.11) 


The results for the first coefficients are: 

C 0 =1; C x = 1; C 2 = 11; 

C 3 = 375; C 4 = 27897; C 5 = 3 817137 

The asymptotic development near tj = » is formulated in the form 

f = f l + f 2 + . . . (9.12) 


where the higher approximations are to be small in comparison with the 
lower approximations, for instance f^ « f^» The first asymptotic 

approximation to correspond to the potential flow is, as was mentioned 
above. 


f 


1 


= TJ - P 


(9.13) 


For this approximation f 1 " = 0 , and one obtains therefore by replacing 
the quantity ff” in the differential equation ( 9 . 8 ) by f 1 f 2 n the 

following equation for the second asymptotic approximation: 
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or 


or 


(tj - 3) t\ + 2f"‘ = 0 



= i O - tj) 
2 


log f" = | Pt) - J T ) 2 + C 


2 

If one substitutes for the integration constant C = — p A + log 7 , 
(7 is thus a new integration constant) one obtains 


f" = 7 e 4 

2 


-f (n-er 


and after integration 


ft 

z 2 


7| e 


- P) 


dr) 


(9.14) 


Because f* («*>)= 1 and f* («) = 0, the solution f* = f* + f* 

1 ' ' 2 v * 12 

satisfies the third boundary condition f* (") = 1. Another integration 

of equation (9.14) gives as second asymptotic solution for f = f^ + fg 


f = TJ 




e ^ dr] 


(9.15) 


This solution still contains two integration constants P and 7 
corresponding to the fact that only one of the three boundary conditions 
was satisfied. The asymptotic solution can, in the same manner, be 
improved still farther by equating f = f^ + fg + f^. The differential 

equation for f^ was solved by Blasius; a more detailed discussion is 

unnecessary . 


These two solutions, the power series near tj = 0 according to 
equation ( 9 . 11 ) and the asymptotic solution near tj = «> according to 
equation ( 9 .I 5 ), now have to be Joined together and the three integration 
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constants a, p and y have thereby to he determined. This is effected 
in the following way: At a point tj = t^, where both solutions are 

serviceable, f, f* and f” from the power series and the asymptotic 
solution are brought to agreement. Because of the differential 
equation ( 9 . 8 ) the higher derivatives will then automatically agree. 

In this manner one obtains three equations for the three unknown 
integration constants. The numerical calculation gives 


a = 0.332; 3 = 1.73; 7 = 0.231 

A table for f, f*, f” taken from a treatise by L. howarth (reference 18) 
is given in table 2*. The velocity distribution in the boundary layer 
u/u o = f'(q) according to equation ( 9 - 5 ) is represented in figure 30 . 

In comparison with the stagnation point profile (fig. l 6 ) it is striking 
that the velocity profile of the plate flow near the surface has a very 
slight curvature. ^At the surface itself it has an inflection point, t hat 

is, for y = 0 : ~ = 0 . 

8 y 2 

The transverse component of the velocity in the boundary layer 
given by equation ( 9 * 6 ) is plotted in figure 31 * It is noteworthy that 
outside of the friction layer, for q-*co 

„ | vn 

T = = 0.865 U c ^ -f 


The fact that on the outer edge of the boundary layer the transverse 
component v / 0 is caused by the deflection of the potential flow from 
the body due to the boundary layer thickness increasing downstream. For 
very large distances from the wall (far in the potential flow) the 
boundary layer solution does not go over exactly into the undistrubed 
potential flow u = U 0 ; v = 0. This has to be tolerated as a (very 
slight) deficiency of the boundary layer solution. 

For the present case a separation of the boundary layer does not 
exist since the pressure gradient equals zero. 


Friction Drag 

From the solution given above the surface friction drag of the 
plate in longitudinal flow is to be calculated. According to equation ( 8 . 9 ) 
the friction drag for one side of the plate is 


* See appendix, chapter XII. 
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t q dx = bpl 


(9.16) 


with h denoting the width, 2 the length of the plate, 
equations (9*7) and (9. 11) 


According to 


(I 11 ) = *l/5 f"(o) = aU 0 \/^2 
\8yA °Vvx 0 V Vx 


Therewith the local surface shearing stress is 


t o - “n u oV^ = 0-3324 u D yS 


(9.17) 


The friction drag according to equation (9.16) is therewith 


W = «|A D 0 \[% J = 


f npiu r 


and therefore the drag of the plate wetted on both sides is 


2W = lwt h \/U 3 nn l 


= 1.328 h\|U 3 up 2 


(9.18) 


If one introduces in the customary manner a dimensionless drag coefficient 
by the equation 


2W 

F f U o 2 


(F = 2 b 2 = wetted area) 


one obtains for the drag coefficient the formula 




(9.19) 
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Displacement Thickness of the Boundary Layer 

By the development of the Boundary layer on the plate, which 
increases downstream with \/x, the potential flow is deflected outward 
from the surface By an amount 5*, which is called the displacement 
thickness of the Boundary layer. It can Be easily calculated from the 
velocity distribution in the Boundary layer, as follows: Let y-j_ denote 

a point outside the Boundary layer; then according to the definition 
for 5* 


or 



( 9 - 20 ) 


According to equation (9*5) 





UT)=0 


[l - f*(Tj)]dTl 




Since the point q = lies outside of the Boundary layer, one can put 

for f ( t) ) the first approximation of the asymptotic solution according 
to equation ( 9 .I 3 ), thus 

f(ni) = Ti x - 3 = m - 1 -'73 


Thus one finds for the displacement thickness of the boundary layer 



6o 
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8 . = 1 . 73 ^* 

o 


( 9 . 21 ) 


The distance from the surface y = 6* is also shown in figure 30. Thus 
the stream lines of the potential flow are, because of the friction 
effect, deflected outward by this amount. 

The actual boundary layer thickness 6 cannot be given accurately 
since the friction effect in the boundary layer thickness decreases 
asymptotically toward the outside. The component of velocity parallel 
to the surface u is asymptotically converted into the velocity U 0 
of the potential flow (the function f’('n) asymptotically approaches the 
value l). If one wants to define the boundary layer thickness as the 
point where the velocity u = O.99 U c (full value), one obtains for it 
according to table 2, q = 5*0. Therewith one has for the thus defined 
boundary layer thickness 



(9.21a) 


The thus defined boundary layer thickness equals about three times the 
displacement thickness of the boundary layer. 

Let here also be introduced the value for the momentum thickness d, 
needed later. This latter is a measure for the momentum loss due to 
friction in the boundary layer and is, as indicated before in equation (6.32), 
defined by the equation 


P" 



J y=o 

The calculation results, because of equation (9.5), i n 



(9.21b) 


Finally the form parameter becomes therewith 
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— = 2.605 
■a 


(9.21c) 


Experimental investigations of the laminar "boundary layer on the flat 
plate were performed by B. G. van der Hegge Zynen (reference 19) and 
M. Hansen (reference 20). In all essential points the theoretical results 
were well confirmed. The measurements showed that the laminar "boundary 
layer exists to about the "Reynolds number (U_ x/v) = 3.5 to 5 X 10^, 

0 crit 

if x denotes the length of run of the "boundary layer. For larger 
Reynolds numbers transition to the turbulent state of flow takes place. 


Fifth lecture on January 5* 19^2. 
"b. The Boundary Layer on the Cylinder (symmetrical case) 

The integration method of Blasius given in the previous section 
was used by Hiemenz (thesis Gottingen 1911) for calculating the boundary 
layer on the circular cylinder. The same method was later further extended 
by Howarth (reference 15) to the general case of a cylindrical body of 
arbitrary cross section. This method will be briefly presented for the 
symmetrical case. One considers (fig. 32) a cylindrical body with 
symmetrical cross section in a flow approaching in the direction of the 
symmetry axis with the velocity U 0 . Let x be the arc length along 

the contour, measured from the front stagnation point, y the vertical 
distance from the surface. Let the potential flow U(x) be given by 
its power series development in x. At the stagnation point (x = 0), 

U(x) = 0, and for the symmetrical case only the odd terms of the power 
series are different from zero. Therefore: 

O R 

U(x) = ux + ux +ux +. . . (9-22) 

13 5 


The coefficients u^, u^, . . . depend Bolely on the shape of the 

body and are therefore quantities known from the potential flow. 


The stationary boundary layer equations according to equation (8.3) 
are also valid for this case with a curved surface and therefore read: 


6u du „ dU , ,, d y u 

dx dy dx dy 2 

|l= 0 

ox oy 




(9.23) 
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From equation (9*22) one obtains for the pressure term: 


U^H = 

dx 


u., 


V 


4u 3 x3 



(9.24) 


The continuity equation is again integrated by the stream function: 


u 



ch|r 

dx 


(9-25) 


It is now necessary to find a suitable formulation for the velocity 
distribution u(x, y), v(x, y) and therewith for the stream 
function i(x, y). In analogy to equation (9.22) a power series in x 
suggests itself for u(x, y) as well, with coefficients, however, which 
are dependent on y. It is important to find a form where the coeffi- 
cients (or functions) dependent on y have a universal character, that 
is, need not be calculated anew for each shape of body, but may be calcu- 
lated once for all. Howarth (reference 15) succeeded in finding such a 
formulation. 


For the distance from the surface one introduces the dimensionless 
variable: 



*(9-26) 


The expression (9.24) for the pressure term suggests that the following 
equation be selected for '4'' 



This yields: (’ = differentiation with respect to tj): 

2 


u = u^xf| + 4u^x^f* + Syp 


u_ 




h 5 


• • • 


(9-28) 


*Qne obtains this equation from that of Blasius according to 
equation (9*3) by substituting for U 0 the first term of the series 

equation (9.22). 
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^ = u f* + 12u x 2 f' + 30 ^^ 
dx ii 3 3” 


u g' + h’ 
5 5 \ 5 


+ • 


(9.29) 


du 

6y 


\ — uif" + 4u x 3 f" + 
I v 1 l 3 3 


6x" 


U 3 

u g" + h" 

5 5 5 


| (9-30) 


>2 u, 
o u _1 

dy 2 ” v 


u i xf> i" + + 6x5 


Vs” + ~ "T 


> ( 9 . 31 ) 


v = — \f— Cu,f, + 12u 0 x 2 f_ + 30 x lf 

hi 11 3 


U 3 

V5 + ^ h 5 


1, 


32) 


After insertion of the expressions (9.24) and ( 9 . 28 ) to (9-32) into the 
first equation ( 9 * 23 ) one obtains by comparison of the coefficients a 
system of ordinary differential equations for the unknown functions 
f^i hej, . . . . , which appears as follows: 


terms with 

u x 
1 

hu^x^ 

c 

6u-. U r I^ 

1 5 

6u 2 x 5 
3 


gives the differential 

f 2 - f f" = 1 + f'" 
111 1 

4f»f» _ 3f"f — f f" = 

13 J 1 3 13 


6 f i g 5 5f l g 5 f i g 5 

6f 'h* - 5f"h - f h" = 

15 15 15 


equation 

1 * f-” 

1 + 8^" 

- + h»" - 8(f» 2 
2 5 3 


f f") 
3 3 


*\ 


> (9.33) 


I 

Formulation of the flow function according to equation (9.27) has 
thus accomplished the elimination of the coefficients depending on body 
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shape (u-p . . . ) from the differential equations- fdr the functions 

f-p fy . . . which thus now have a universal character. 

The houndary conditions for the functions f^, f^, . . . follow 

from 


y = 0: u=v=0 

y = u = U 


hy comparison of equation ( 9 . 28 ) with equation ( 9 . 22 ); they are 



h 


5 



= 0 ; 



( 9 - 34 ) 


The differential equations (9*33) are all of the third order, and 
equation (9.3I4.) gives tlrcee “boundary conditions for each. The differen- 
tial equation for f,(T|) is non-linear and is identical with the differ- 
ential equation (6.28) obtained in chapter YI for the stagnation point 
flow: f-^ s qpj r\ = | , as follows by comparison of equation ( 6 . 26 a) 

with equation ( 6 . 26 ). All the remaining differential equations in 
equation ( 9 . 33 ) are linear, with the coefficients determined by the 
functions of the preceding approximations. 

The solutions of the differential equations (9*33) are best obtained 
by numerical integration. The functions f^ and f^ were already 

calculated by Hiemenz (reference 10). Howarth (reference 15) improved 
the tables for f^ and recently Nils Frossling (reference l 6 ) calculated 
g and h as well. The f! which is identical with cp 1 according to 

5 5 1 

equation (6.28) was already given in figure l 6 . The function f^ can be 
seen from figure 33* 'the functions g^ and h^ from figure 3^» The 
numerical values are compiled in table 3 - 

Concerning the applicability of this calculation method it must be 
mentioned that for slender body shapes the series for U(x) and u(x,y) 
converge poorly. The reason is, that for these body shapes U(x) has 
a very steep ascent in the neighborhood of the stagnation point (fig. 35 )* 
while showing a rather flat curve further on. Such a function cannot be 
developed readily into a Taylor series. For such body shapes many more of 
the functions of the differential-equation system (9*33) would be required 
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than have "been calculated bo far. For hlunt bodies, as for instance the 
circular cylinder, the convergence is considerably better so that one 
proceeds rather far with this calculation although not always up to' the 
point of separation. 

Howarth (reference 15) also performed the corresponding calculation 
for the unsymmetrical case where the even coefficients also appear in 
the power series for U(x). This is the case for a symmetrical body 
at an angle of attack and quite generally for any unsymmetrical body. 

Frossling (reference 16) made the application to the rotationally 
symmetrical case. 


Circular Cylinder 


The boundary layer on the circular cylinder will be treated as an 
example for the application of this method. Whereas Hiemenz (reference 10) 
took a pressure distribution measured by him as the basis for this case, 
we shall here calculate with the potential— theoretical pressure distribution. 
The velocity distribution of the potential flow reads, with the symbols 
according to figure 3 6, 


U(x) = 2U 0 sin <p = 2U 0 sin f (9-35) 

X\ 


The power series development gives: 


U(x) = 2 U. 


£ - -L (Z ) 3 + -L + 

R 3! \R/ 51 VR / 


(9.36) 


In comparison with equation (9.22): 

U U TJ ' 

u, = 2 — ; u, = — — u,- = — — ; 

1 R 3 3* 5 5! ~5 


55 R 5 


(9.36a) 


Therewith follows from equation (9.26): 



(9-37) 


It follows that for the velocity distribution from equation (9.28) 
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1 u(x,y) = x , 

2 U 0 R 1 



f* 

3 


+ 20h.)- 


(9-38) 


One further calculates the position of the 
is, according to equation (8.5) determined 

results from equation (9*38): 


separation point 


by 


du 

8y 


= 0 . 


y=o 


x^ which 
Therewith 


f"(0) ^ f"(o) 

1 R 3! 3 



. . = 0 


(9.38a) 


With the numerical values f^(0) = 1. 2326b; 


f^'(0j = 0.72b6 


one finds: 


X 

— = 1.60; <p = 92° 

R A 


(9-39) 


Hiemenz (reference 10) "based his calculation on his experimental 
pressure distribution; he calculated the separation point to be at 

cp = 82°, whereas his measurements gave cp = 8l°. This result is 
A A 

considerably different from that obtained with the potential— theoretical 
pressure distribution. The reason is that for a body as blunt as the 
circular cylinder the experimental and the potential— theoretical pressure 
distribution in the neighborhood of the separation point differ greatly. 

The method described here of calculating the boundary layer by a 
power— series development starting from the stagnation point has found 
but little acceptance because of the extensive calculation required. 

It is, however, indispensable for fundamental considerations, since there 
exist no other exact solutions of the differential equations of the 
boundary layer for the flow about a body. 

Thus approximation methods came into use for the practical per- 
formance of boundary layer calculations; they will be discussed in the 
following chapter. It is true that their accuracy is sometimes consider- 
ably lower than that of the previously discussed exact solutions. 


c. Wale behind the Flat Plate in Longitudinal Flow 

The application of the boundary layer equations is not absolutely 
limited to the presence of solid walls; They may also be applied if 
there is present within the flow a layer in which the friction effect is 


*This result varies somewhat if in the series development 
equation (9.38a), one tales further terms into consideration. However, 
for this purpose one would have to calculate at least up to the term x 
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predominant. This is the case for instance when within the flow two 
layers of different velocities adjoin, as for instance in the wake region 
behind a body or at the outflow from an opening. In this and the following 
chapter we shall treat two examples of such flows which we shall later 
encounter again in the discussion of turbulent flows. 

4 

The wake flow behind the flat plate in longitudinal flow is chosen 
as the first example (fig. 37)* At the trailing edge of the plate the 
two boundary layer profiles grow together and form a "wake profile" the 
width of which increases wAh distance while the velocity decrement at 
its center decreases. The size of the "wake" is directly connected with 
the drag of the body. Otherwise, however, the form of the velocity 
distribution in the wake at a large distance from the body 1 b not dependent 
on the shape of the body, whereas the velocity distribution very close 
behind the body naturally depends on the boundary layer of the body and 
on any existing separation. 

From the velocity distribution measured in the wake one may cal- 
culate the drag by means of the momentum theorem in the following manner. 
The momentum theorem states: The variation of the momentum with time 

( = momentum flow through a control area fixed in space) equals the 
sum of the resultant forces. By resultant forces one has to understand: 

1. Pressure forces on the control area, 2. Extraneous forces, which are 
transferred from solid bodieB to the flowing fluid, for instance the 
shearing stress at the surface which gives the friction drag. For the 
present case the control area AAjBB-^ is placed as indicated in figure 38 . 

Let the boundary A^B^ which is parallel to the plate be so far distant 

from it that it lies everywhere in the undisturbed velocity U 0 . Further- 
more, the rear' cross section BB-j_ 1b to lie so far behind the plate that 

the static pressure there has the same undisturbed value as in front of 
the plate. Then the pressure is constant on the entire control area, so 
that the pressure forces make no contribution. In calculating the 
moment 'am flow through this control area one has to consider that, due to 
continuity, fluid must flow out through the boundary A-^B^, namely the 

difference between the larger quantity flowing through the cross 
section AA^ and the smaller quantity flowing through the cross 

section BB^. The cross section AB does not make a contribution to the 
x-momentum, since for reasons of symmetry the transverse velocity on it 
equals zero. The momentum balance is given in the table below, with 
entering momentums counted as positive, outgoing ones as negative. 
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Cross section 

1 

Mass 
m3 /sec J 

.. . — - - 

Momentum flow in x— direction 
= volume x density x velocity 

A B 

0 

0 

A A^ 

b 

«c 

ih 

u 0 dy 

• p b 
1C 

|h 

Uo 2 ^ 

B B x 

- b 

u: 

Th 

u dy 

- p b 

J< 

ih 

u 2 dy 

A 1 B 1 

- 1 

i 

nh 

^ (U 0 - u)dy 

Jo 

- P b U ( 

nh 

j ,(U C - u) dy 

Jo 

1= 

Control area 

Mass 
= 0 

y~ Momentum Flow 
= W 


The total momentum loss of the flow for the present case equals the 
drag W of -one side of the plate. Thus one obtains 


/ft 


n°° 




u (U_ - u) dy = b J 


Uo 


u (U Q - u) dy 


UO 


(9.40) 


The integration therein may be extended from y = 0 to instead of 

from y = 0 to h, since for y > h the integrand in equation (9.40) 
vanishes. For the drag of the plate wetted on both sides one obtains 
therefore : 
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l +eo 


2 W = bp 


u(U Q -u) dy 


(9.41) 


In equations (9*40) and (9*41) the integrals are to "be extended over the 
wake as indicated at a distance "behind the plate where the static pressure 
has its undisturbed value. However, one may naturally apply equations 
(9.40) and (9.41) also to the boundary layer on the plate at a certain 
point x; then they give the drag of the part of the plate from the 
leading edge to this point. The definite integrals in equations (9*40) 
and (9*41) represent physically the momentum loss due to the friction 
effect. As mentioned before, it is customary to introduce for this 
integral also the momentum Iobb thickness d by the following equation 
(compare equation ( 6 . 32 )). 

U Q 2 18 = [ u ( u 0 - Ti) dy (9.42) 


7=o 


Therewith the formula for the drag may also be written, by comparison 
with equation (9.40): 


W = b 


2 


P u 0 0 


(9.43) 


The velocity distribution in the wake, particularly at large 
distance x behind the plate in longitudinal flow (fig. 37 ) is to be 
calculated next. ThiB calculation must be performed in two steps: 

1. By a development "from the front", that is, by a calculation which 
follows the further development of the Blasius-boundary layer profile 
present at the trailing edge of the plate. 2. By a development "from 
the rear" that is, by an asymptotic calculation for the wake, under the 
assumption that the difference velocity 


u-^x, y) » T7 0 - u(x, y) 


(9.44) 


i s small . 

The first calculation was performed by S. Goldstein (reference 21) 
and is very troublesome; the second'was indicated by Tollmien (reference 3 ) 
and yields rather simple laws, in particular also an exact solution of the 
differential boundary layer equation. Since such exact solutions are 
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comparatively rare and since, moreover, the asymptotic law for the wake 
applies not only for the flat plate in longitudinal flow hut for any 
arbitrary body shape, this asymptotic solution will be treated here' 
somewhat more thoroughly. The wake velocity u^(i, y) introduced in 

equation (9.44) is assumed to be so small in comparison with the free 
stream Telocity U 0 that the second— power terms (u-^/Uq) 2 are negligible 

relative to 1. Moreover, the pressure term dp/dx in the boundary layer 
equation is to be set equal to zero for the first asymptotic approximation. 
Therewith the differential boundary layer equation (8.3) becomes, for the 
present case: 


XL 




(9.45) 


With the boundary conditions: 


du-i 

y = 0: - — = 0 

oy 

y = »: VL ± = 0 (u = U o ) 


(9.45a) 


For the solution of the differential equation (9.45) one introduces as 
before in the plate flow according to Blasius* equation (9*3) the new 
variable 


n 



(9.46) 


Further, one uses for 


u 


the equation: 


_1 

U 1 = u o c (!) 2 s(ti) 


(9.47) 


The distance x from the trailing edge of the plate is thus made 
dimensionless by dividing by the plate length. The power — ^ for x 

ia given by the fact that the momentum integral which, according to 
equation (9.41) gives the plate drag must be independent of x. With 
the second— power terms neglected the drag of the plate Wetted on both 
sides is, according to equation ( 9 . 41): 
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ft + oo 


2V = b p U D 


\ dy 


(9.48) 


y=— oo 


One finds further: 


1+oi 


^ + 00 


n+» 


V X ' y)dy = U Q C (|) M g(t|)dij = 


u c \fil 
° V U o 


g(Ti)drj (9 .49) 


and therewith: 


2W = h p 


u Q 2 cySl g (l) ) 4, 


(9-50) 


The calculation of the single terms in equation (9.45) gives: 

\ 

sr = ~ u o I (f) 2 | (g + T) g ») 


8u, 


^ = u ° c (f) 2 s ' 


O U 1 /yN 

v TV - 0 (!) 

oy 


1 

2 1 


x 8 




(9.51) 


By insertion into equation (9*45) one obtains after division by 

2 — — i 

0 Uo ( X A) ^ x the following differential equation for the velocity 
distribution g(Tj): 


g" + | Tig* + | g = 0 


(9-52) 
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with the boundary conditions: 

T) = 0: g* = 0; T| = oo : g = 0 

A single integration gives: 

g» + i T| g = K = 0 


where, because of the boundary condition at q « 0, the integration 
constant K must be zero. Repeated integration gives the solution: 



(9-53) 


where a multiplicative integration constant may be put equal to one 
without limiting the generality since the velocity distribution Uj 

still contains, according to equation (9.47), the multiplicative free 
constant C. This constant C is determined from the consideration 
that the plate drag obtained from the momentum loss (equation (9*50)) 
must be the same as the frictional drag of the plate. First, 


n+oo 


r +0# _ n_ 


g(Tl) dT) = 


(J-e 


e dq = 2 \J~a 


and therewith from, equation (9.50): 


2W = p b U 
H o 


2 



On the other hand, the friction drag of the plate wetted on both sides 
was according to the solution of Blasius (equation (9.18)) 


2W = I.328 




Therefrom follows 2C \[n = I.328 and 
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0.664 

\T*~ 


(9.5*0 


Thus the final solution for the wake velocity for the flat plate in 
longitudinal flow "becomes: 


_ i i y 2 u c 

1 _ 0 . 664 / x \ 4 x v 

u 0 -\pr [ij 


(9-55) 


The velocity distribution of this asymptotic law is represented in 
figure 39. It is noteworthy that the function for the velocity distribution 
is identical with Gauss* error function. In keeping with the hypothesis 
the law according to equation (9*55) is valid only for large distances 
behind the plate, according to the calculations by Tollmien (reference 3) 
for x - 31. 

The development of the wake from the front, performed by Goldstein 
is valid only for comparatively small x/l . However, for intermediate 
x/l both solutions can be joined to some extent, so that one then obtains 
the velocity distribution in the entire wake. Such a figure is given by 
Tollmien (reference 3). 


Sixth Lecture on January 12, 19^2. 
d. The Plane Jet 

A further example of a flow without bounding wall to which the 
boundary layer theory is applicable is the outflow of a jet from a hole. 
The problem to be treated is the plane stationary one where the jet goes 
out from a long narrow slot and mixes with the surrounding fluid at rest. 
This is one of the rare cases where the differential boundary layer 
equations may be integrated exactly. The calculations were performed by 
H. Schlichting (reference 22) and W. Bickley (reference 30) and will be 
discussed a little more thoroughly. 

Due to the friction effect the jet entrains a part of the fluid at 
rest and sweeps it along. There results a stream— line pattern like the 
one drawn in figure ^0. Let the x— direction coincide with the jet axis 
and the origin lie in the slot. It then immediately becomes clear that 
the width of the jet increases with the distance x and the mid— velocity 
decreases with the distance x. For the calculation the slot is assumed 
to be infinitely narrow. In order to make the volume of flow, together 
with its momentum, finite, the velocity in the slot is then infinitely 
large. Again, as in the previous example concerning the wake flow, the 
pressure term dp/dx may be neglected since the pressure varies only 
very little in the x— direction. The smallness of the pressure term can 
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also be shown subsequently, from the finished solution. Thus the differ- 
ential boundary layer equations for the present case read, according to 
equation (8.3): 





du 

cbc 



0 


*\ 


f (9.56) 


with the boundary conditions: 


y = 0: 



y = 00 : u = 0 


> 


* 


(9.56a) 


Since the pressure is constant, the entire momentum flow in the x— 
direction for control area fixed in space (compare figure 40) must be 
independent of the distance from the hole x. If one chooses the lateral 
boundaries of the control area at so large a distance from the axis that 
there u = 0, then 


p +00 


J = p 


u 2 dy 


= Constant 


(9.57) 


It is to be noted for the integration of the equation of motion 
equation (9.5 6) that for this problem, as before for the plate in longi- 
tudinal flow, no special length-dimension exists. Thus affinity of the 
velocity profiles u(x, y) is suggested, that is: with b signifying 

a suitable width of the Jet, the velocity profiles are only functions 
of y/b. Accordingly one uses the following expression for the stream 
function \|r: 




(9.58) 
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The two - at first unknown - exponents p and q are determined from 
the conditions that 

1. the momentum flow for the x-direction is independent of x 

according to equation (9*57), and that 

2. the acceleration terms, for instance u — , and the inertia 

dx 

term in equation ( 9 * 56 ) are of the same order of magnitude 
and hence must he of the same degree in x. 

This yields the two equations 


It follows that 


2p — q - 0 
2p-2q-l = p-3q 


> 


V 


(9-59) 


p = 1/3 } q = 2/3 (9.60) 


Therewith the final equations, after addition of suitable constant factors, 
read as follows: 


T) = — — 

1/2 2/3 
3v i 


(9.61) 


♦ = V l/2 x l/3 f ( tj ) (9.62) 

Therefrom one obtains, with 

= _ £ a. M 

dx 3i’ Sy l/2 

3v 


the following expressions for the velocity components and their derivatives 
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u 


v 




3* 


1/3 


f'(Tl) 


- i v 1//2 I 2 / 3 (f-2nf») I 


( 9 . 63 ) 


s 


N J ‘/3 
S = - X 9 

- I_1 | v " 1 / 2 f” 


V 



-5/3 


_ ft *» 

27 


By substitution into the differential equation ( 9 . 56 ) there results, 

1 -5/3 

after cancelling the factor — x , the following differential equation 
for the flow function f(rj): 

f,2 + ff n + f .t» - o (9.64) 

with the boundary conditions: 


rj = 0 : f = 0 ; f" = 0 

tj = 00 : f» = 0 


(9.64a) 


As for Blasius* plate flow here also an ordinary differential equation ( 9 . 64 ) 
was obtained from the two partial differential equations ( 9 . 56 ) by means of 
the similarity transformation equation ( 9 . 61 ). Here also, as in most 
boundary— layer problems, the differential equation is non-linear and of 
the third order. 

The integration of this differential equation (9.64) is attained in 
a surprisingly simple manner. First, one obtains by a single integration 


ff* + f" = Constant = 0 


(9.65) 
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The constant of integration therein ie zero because of the boundary- 
condition f"(o) = 0. The second order differential equation now obtained 
could be integrated once more if a factor 2 were present in the first 
term. This factor can be obtained by performing the following further 
similarity transformation: 

One puts: I = a t) (9.66) 

f - 2a F(6) (9.67) 


a is a free constant which is determined later. With the equations (9.66) 
and ( 9 . 67 ) one obtains from equation ( 9 . 65 ), the prime now signifying 
differentiation with respect to 


F" + 2FF» = 0 


( 9 . 68 ) 


with the boundary conditions: 


1 = 0: F = 0; £ = »: F* = 0 


(9.68a) 


This differential equation can now be integrated again and yields 


F* + F = K 


The constant of integration K is obtained from the boundary conditions, 
equation ( 9 . 68 a), as K = 1, if one stipulates F'(o) = 1, which is 
possible without limiting the generality because a is still present as 
a free constant in F. One now has for F the first order non-linear 
differential equation 


F 2 = 


(9.69) 


which is a Riccati differential equation. The integration yields 


I = 



1 1 + F 

1 log = arc tanh F 

2 1 — F 


and therewith for the inverse function 
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-2t 

F = tanh I = i-Z_® — - (9-70) 

-25 

1 + e 


Furthermore, ~ = 1 — tan^ | . if one inserts the solution found into 
dt 

equation ( 9 * 63 ), one obtains for the velocity distribution 


_ 1 

u = j x — tanh^l) 


(9.71) 


The velocity distribution over the width of the jet calculated from this 
equation is represented in figure 4l. The free constant a remains to 
be determined. This can be done from the condition (equation (9*57)) 
according to which the momentum flow in the x— direction is independent 
of x. From equation (9.71) and (9.57) one obtains 


r 


J = | p a3v X / 2 (l - tanh 2 ! ) d| = y p c^v 1 / 2 

Jo 


(9.72) 


Let the momentum J for the Jet be a prescribed constant which is, for 
instance, directly related to the excess pressure under which the jet 
flows from the slot. Then a becomes, according to equation (9.72), 


a = 0.8255 


j y/3 

172) 

pv / 


(9.73) 


and therewith the velocity distribution 


U = 0.4543 I 

(y 

1 2 

r. 

f 1 — tanh 2 1 ) 
\ 


\p vx 

J 



Jv ^ 

2 , 

1/3 

j » ^ V 

V = 0.5503 ( 

) ’ 

2£ (1- tanh^l ) - tanh 1 


VP* J 


l J 



> 





(9.74) 


I = 0.2752 
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The value of the transverse component of the velocity at the edge of the 
jet (y = <») also is noteworthy. From equation (9.7h) one finds for 
this lateral inflow velocity 




= - 0.550 




( 9 - 75 ) 


One can further calculate the flow volume for a layer of unit height 



One finds 


£ 

P 



(9.76) 


The fltiw volume increases downstream, since fluid at rest is carried 
along from the side. 


The solution indicated here naturally has a singular point at 
x = 0, since an infinitely narrow slot with infinitely large exit 
velocity was assumed. Actually for a narrow slot of finite width one 
has immediately behind the slot opening a velocity distribution that is 
rectangular across the jet cross section but which at some distance is 
transformed into the bell— shaped distribution found here with width 


b x 


2/3 


and mid— velocity ~ x 


~ 1/3 


Finally it should be mentioned that the corresponding rotationally— 
symmetrical problem where the Jet goes out from a very small circular 
hole also can be solved in closed form, (compare H. Schlichting 
(reference 22)). In this case the width of the jet is proportional to x 
and the midvelocity proportional to x~ . 


e. The Boundary Layer for the Potential Flow U = u^x m . 

Another class of exact solutions of the boundary layer equations 
will be discussed briefly which includes the plate in longitudinal flow 
and the stagnation point flow as special cases. Falkner and Stan 
(reference 37 ) have shown that, just as for Blasius* plate flow, the 
boundary layer differential equations for the potential flow 


U(x) = u^ x m 


( 9 . 77 ) 


80 


NACA TM No. 1217 


can "be reduced by a similarity transformation to an ordinary differential 
equation (u^ = constant, m> 0 accelerated, m< 0 retarded flow). For 
m > 0, x = 0 is the stagnation point of the flow. For m = 0 one 
obtains U = u^ = U Q , therefore the plate flow; m = 1 gives U = u^Xj 

therefore the stagnation point flow according to equation (6.17). 

The differential equations of the boundary layer read 

8u . „ 8u TT dTJ d 2 u 

u — + v — = U — + v 

dy dx 3y 2 

** + = 0 
8x 8y 


The pressure term becomes 


U dll = m u 2 

dx 1 


2m-l 

x 


As a new independent variable one introduces 

mr-1 



(9-78) 


and the continuity 
function for which 


equation is integrated by introduction of a stream 
one uses the equation: 



(9-79) 


One has 


81 _ m — 1 
dx 2 X* 


m-1 



and one obtains 
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m 

u = qp'(l) = Wp»(l ) 


( 9 . 80 ) 


du 

dy 



n q-m-g 

bn m + 1 p ~p~ . 

V - — = u , x d cp«"(|) 

2 2 1 Y V5 ' 



r m— 1 


cp* + 111 — - 


I 



v 



1 



m + 1 


<P + 




( 9 - 81 ) 


After insertion into the first equation of motion and division by 
m u 1 2 x 2n ^ 1 one obtains, when 


2m _ 
m + 1 


( 9 . 82 ) 


the following differential equation for <p(|): 


cp*" = - qxp" + 0 



( 9 - 83 ) 


Boundary conditions: 


|=0: cp = qp = 0; | = 00 : qp* = 1 

The differential equation (9.83) was solved for different values by 
Hartree (reference 38). The result is given in figure 4 la. The 
corresponding values of 0 and m are given in the following table. 
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in 

p 

-0.0654 

- 0.14 

0 

0 

0.111 

0.2 

0.333 

0.5 

1 

1 

k 

1.6 


For accelerated flow (m > 0 , B > 0 ) one obtains velocity profiles 
without inflection points, for retarded flow (m< 0 , 3 < 0 ) velocity 
profiles with inflection points. Separation occurs for 


3 - - 0.199, that is, m = - 0.091 


Separation takes place for the potential flow U(x) = u thus 

for very weak retardation. Compare chapter XI a, where^almost the same 
result is obtained from an approximation calculation. 


Special cases: 


1 . Stagnation point flow: It is obtained for 


m = 1; 3 = 1 


Then | = 





— y; 



x q>(6). 


These are the same expressions as 


for the stagnation point flow, equation (6.26a) and ( 6 . 27 a), also (6.19) 
and ( 6 . 26 b). The differential equation (9.83) also is transformed into 
the equation of the stagnation point flow (equation (6.28)). 


2 . Plate in longitudinal flow: This case is obtained for 


m = 0 ; 3=0 
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Then 



7 = — \| — y = 4 i with 


\fi 


q signifying Blasius* 


variable according to equation ( 9 . 3 ). Furthermore, + becomes 

compared with the expression 


t = \[? \/v U qX cp(£); thus <p = ~ 

for \|r for the plate flow equation ( 9 . 4 ). Because of — = — , 

r~ ^ dT l 

<p"(£) = f M (q), and <p ,w (S) = 2 f'"(q), the differential equation ( 9 . 83 ) 

is for this case transformed into 2f , "(T|) + ff”(q) = 0. This is identical 
with equation ( 9 . 8 ). 


CHAPTER X. APPROXIMATE SOLUTION OF THE BOUNDARY LAYER EQUATION 
BY MEANS OF THE MOMENTUM THEOREM 
(KARMAN-POHLHAUSEN METHOD, PLANE PROBLEM) 
a. The Flat Plate in Longitudinal Flow 


The examples of exact solutions of the boundary layer equation 
discussed in the previous chapter give sufficient proof of the rather 
considerable mathematical difficulties in solving the differential 
equation. Yet the examples treated were selected as simple as possible. 

In some other cases the mathematical calculations are still more difficult 
than in those, examples. Particularly the problem of flow about a body 
of arbitrarily prescribed shape cannot, in general, be solved by exact 
solution of the differential equations of the boundary layer. Just this 
problem, however, is of special practical importance, for instance for 
the calculation of the boundary layers on airplane wings. 

There exists therefore a tendency to apply at least approximate 
methods, even if their accuracy is sometimes not quite satisfactory, 
for cases where the exact solution cannot be obtained with a reasonable 
expenditure of calculation time. Such simpler approximate solutions can 
be attained if one does not attempt to satisfy the differential boundary 
layer equation for every particle of fluid. Instead one selects a 
plausible expression for the velocity distribution in the boundary layer 
and satisfies merely the momentum equation which is obtained by inte- 
gration from the equation of motion. A method on this basis for the 
plane problem- of flow about an arbitrary body was worked out by von Karman 
and Fohlhausen (references 23 and 24). We shall try out this method in 
this chapter at first on the simpler case of the flat plate in longi- 
tudinal flow, where no pressure variations exist. For this special case 
the momentum theorem yields the statement that the momentum loss of the 
flow through a control area fixed in space according to figure 42 equals 
the friction drag W(x) of the plate from the leading edge (x = 0) to 
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the point x. Application of the momentum theorem for this case was 
discussed in detail in chapter IX; for the drag of the plate wetted on 
one side according to equation (9.40) it had resulted in the formula: 


nee 


W(x) -Dp 


u (U Q - u) dy 


( 10 . 1 ) 


dy=0 


On the other hand the friction drag can also he expressed as the integral 
of the surface shearing stress, namely: 


W(x) = h t d (x) dx 


ix=0 


( 10 . 2 ) 


In forming the integral (equation (10.1)) it is to be noted that the 
integrand outside of the boundary layer, where u = U Q , does not make 

a contribution. By comparison of equations (10.1) and (10.2) it follows 
that: 


nee 


T 

O 



u (0 o 


Uy=o 


u) dy 


(10.3) 


If one introduces moreover the momentum loss thickness, as defined in 
equation (9.42), equation (10. 3) can also be written in the form: 


tt 2 _ lo 

° dx " 0 


(10.4) 


This is the momentum theorem of the boundary layer for the special case 
of the flat plate in longitudinal flow. Physically it states that the 
surface shearing stress equals the momentum loss in the friction layer, 
since in the present case the pressure gradient makes no contribution. 
The next chapter will acquaint us with the extension of equation (10.4) 
to include the more general case of a boundary layer tfith pressure 
difference. 

Equations (10.4) and (10.3), respectively, will now be used for 
approximate calculation of the friction layer on the flat plate in 
longitudinal flow. A comparison of the results of this approximate 
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calculation with the exact solution according to chapter IXa will give 
information about the usefulness of the approximation method. For the 
approximate calculation one selects a suitable expression for the 
velocity distribution in the boundary layer in the form: 

u-U 0 f(l)- u 0 f(i) (10.5) 

with 

t) = 5 = 5(x) (io.6) 

o 


5 represents the boundary layer thickness, undetermined at first. For 
the flat plate it may, moreover, he assumed again that the velocity 
profiles at various distances from the leading edge of the plate are 
affine to each other. This assumption is contained in equation (10.5) 
if f(n) there stands for a function which no longer contains any free 
parameters. Furthermore, f(T]) must, for large values t), assume the 
cons t ant value 1 . 


The velocity distribution being given by equation (10. 5)^ the 
momentum integral in equation (10. 3) may be evaluated. It yields: 


u (n - u) dy = u 


f (1 - f) dn 


1=0 


(10.7) 


The definite integral in equation (10.7) can be calculated immediately 
if a definite formulation is given for f ( T) ) . Thus one puts 


a = 


f (1 ~ t) <1t) 


UT|=o 


( 10 . 8 ) 


Hence 


n~ 


2 

n o «- 


u(U — u) dy = a 6 U 


(10.9) 


Uy=o 


or: 


d = a 6 


( 10 . 10 ) 
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Furthermore, the displacement thickness of the boundary layer 6* becomes, 
according to equation (9.20): 


n 1 


6 * = 6 


(1 - 


f) dn = 0 ,^ 5 


(10.10a) 


lio 


On the other hand, the shearing stress t q at the surface is: 


1° = V f * (0) = 3 ^2- 


dy/ 8 

y=0 


( 10 . 11 ) 


if one introduces the further simplification 


3 = f‘(0) 


( 10 . 12 ) 


By introduction of these expressions into the momentum equation (10.4) 
there results: 


or 




(10.13) 


5 d5 = il 

dx a U 

o 


(10.14) 


The integration with the initial value 6=0 for x = 0 yields, as 
first result of the calculation: 


6 


\/Tf Jvx 

V « V U o 


(10.15) 


Hence the shearing stress becomes, according to equation (10.11): 



(10.16) 
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and furthermore 

\ l 

T o dx = 

Jo 


and hence, finally, the total drag of the plate wetted on both sides 
according to equation (10.2): 



(10.17) 


A comparison of the results for boundary layer thickness, surface shearing 
stress, and total drag, which were thus found, with the corresponding 
formulas for the exact solution according to equation (9« • •) shows that 
the approximate calculation according to the momentum theorem reproduces 
the characteristics of the formulas with perfect correctness in all cases, 
that is, the dependence of the boundary layer quantities on the length 
of run x, the free stream velocity U 0 , and the viscosity coef- 
ficient V. The numbers a, p, unknown at first, can be determined only 
after making special assumptions for the velocity distribution, that is, 
explicitly prescribing the function f(t)) in equation (10.5). 


Numerical example s 


The usefulness of the method of approximation will be investigated 
by a few numerical examples. The accuracy of the results will depend to 
a great extent on a suitable choice of the expression for the velocity 
distribution according to equation (10. 5). At any rate the function f(q) 
must equal zero for tj = 0 and have the constant value 1 for large t| . 
As first example we select the very rough assumption that the velocity 
distribution in the boundary layer is represented by a linear f unction 
according to figure l*3a. 


Thus : 


0 — i) — 1: f(n) = 7 

tj - 1: f(n) = 1 

( 10 . 18 ) 


Hence the results for the numbers 
and (10.12) are: 


P according to equation (10.8) 



P = 1 


(10.19) 
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The formulas (10.15)* (10.16) and (10.17) can now he evaluated immediately. 
One obtains the results: 


8 = 2 \J3 = 3-464 ( 10 . 20 ) 

T o = -3j= uU 0 \[^= 0.289 (10.21) 

2 V Vx v Vx 

2W = h upU 0 3 X = 1.155 b \| npU 0 3x (10.22) 


A velocity distribution in the form of a cubic parabola according to 
figure 13b is selected as second numerical example in the following manner : 


0 < n < 1: f(n) 

t) - 1: f(Ti) = 1 


(10.23) 


This satisfies the conditions: 


T) = 0: 


f = 0; 


T| = 1: f = 1; 


f» = 0; 


that iSj ths boundary layer profile joins the velocity of the potential 
flow with a continuous tangent. The calculation of the numerical factors 
according to equations (10.8) and (10.12) gives: 

a - 2$; » = I < 10 - 21( > 


and hence for the characteristic parameters of the boundary layer: 


5 = k.6b\[V£ 



(10.25a) 


T o = °* 323 uU 0 


(10.25b) 
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2W = 1.29 1 \J hpUq 2 


The exact value for the drag is, according to equation (9.18), 

2W = 1.328 "b ^U 0 VpI. The simple assumption of a linear velocity distri- 
bution therefore gives a drag too small hy thirteen percent, the cubic 
velocity distribution a drag too small by three percent. 

The calculation of the displacement thicknesses of the boundary 
layers according to equation (10. 10a) results, for the linear velocity 

distribution, in 5* = i 5, and for the cubic velocity distribution 
3 ^ 

in 5* = g 5. This gives, because of equations (10. 15) and (10.25a): 


6* = 1.732 


(linear velocity distribution) 


6* = 1.7^0 (cubic velocity distribution) 


(10.26) 


The agreement with the exact value 6* = 1.728 y— according to 

equation (9.21) is surprisingly good; this is, however, more or less 
accidental. 

The essential characteristics of the boundary layer according to 
the approximate calculation described above are once more compiled with 
the exact solution in the table below. 


Characteristics of the Boundary Layer on the Flat Plate; 
Comparison of Approximate Calculation and Exact Solution 


Kind of calculation 5 



Linear approximation 1.732 0.578 

(fig. k3a) 

Cubic approximation 1.7^0 0.6k 5 

(fig. h3b) 


Exact solution 
(Blasius) 



-729 
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As one can see from this table, the agreement, particularly of the cubic 
approximation and the exact solution, is rather good. On the whole, the 
results of this calculation with the aid of the momentum theorem may be 
regarded as very satisfactory, especially in view of the simplicity, as 
compared to the exact calculation. 


Seventh Lecture (January 19, 1942) 

b. The Momentum Theorem for the Boundary Layer 
with Pressure Drop (Plane Problem) 

Last time the boundary layer on the flat plate in longitudinal flow 
was calculated by means of the momentum theorem. Today the general case 
of the boundary layer with a pressure difference' in the flow direction 
will be treated. One considers the flow along a curved surface, and 
measures the coordinate x as arc length along the surface; let y be 
the perpendicular distance from the surface, U(x) the given potential 
flow (fig. 44) . The fundamental equation may be obtained by a momentum 
consideration as in chapter IXc; now, however, the pressure difference 
has to be taken into consideration. The same result is obtained by 
integration of the equation of motion of the boundary layer with respect 
to y from y = 0 (surface) to y = h, the layer y = h lying 
everywhere outside of the friction layer (fig. 44). 


The differential equations of the boundary layer for the steady 
case read, according to equation (8.3)* 


du | 8u 1 dp 

dx by p bx 

5u , dv _ q 
bx by 



(10.27) 


with the boundary conditions: y = 0: u = v = 0; y=»: u = U. The 

integration from y = 0 to h gives: 


n h 


Oh 


1 _d 

2 dx 


^ dy + 


by dy 


hi^+V 
p dx 


du 

dy 


(10.28) 


Lr J JO 


Jy=0 


In the first term the differentiation with respect to x and the integra- 
tion with respect to y may he interchanged, since the upper limit h 
is independent of x. On the left side the second term is transformed by 
integration by parts: 
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f»h 


T | dJ 


' l h 
u v 

L JO 


Oh 


^udy = U 

py 


Vu - 


t udj 


Uo 


T h re P re renting the transverse velocity outside of the boundary layer. 

By continuity, — = — ^ and 
dy 


Oh 


\ 




( 10 . 29 ) 


and hence: 


Oh 


_ du . 
v — dy = 

dy 


- u 


Uo 


du . . 


UO 


du 

dx 


ay 


Uo 


Insertion in equation (10.28) gives, "because of: 

ih 


fv T 

du _ o 

Mo p 


the relation: 


Oh 


dx 


dy - U 


ay = _hd£_Io 


dx 


p dx p 


Jo 


(10.30) 


(10.31) 


/ / 

This is the so-called Karmen integral-condition , first given by v.'Karman 
(reference 23). 

For the pressure term one now introduces the potential velocity U(x); 
furthermore, equation (10.31) is to be transformed so that the displacement 
thickness 5 * and the momentum thickness -3 appear in it as defined by 
equation (6 . 31 ) and (6.32), namely: 
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U S* = (U — u)dy 


U •a = u (U - u) dy 


According to Bernoulli’s equation: 


1 dp = _ u au 

P dx dx 


dU dU 

+ U — 

dx dx 


which can also "be written: 


= U 2 dy + U U dy 

p dx dx dx J 


By substitution of equation (10. 3*0 into equation (10. 31) there results: 


— = — (U 2 - u 2 ) dy - U — (U - u) dy 

p dx dx 


(U-u) u dy + (U — u) dy / — U •— (U - u) dy 


and after differentiation of the second term; 


— = — (U — u) u dy + — (U — u) dy 
p dx dx 
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The displacement thickness and the momentum thickness can now be introduced 
directly and one obtains: 

— = — (U 2 ^) + 32 6* u 

p dx dx 


or 


(10.36) 
J 




6*) U SI 

' dx 


This is the form of the momentum equation for the boundary layer with 
pressure drop that will be used as a basis for further considerations. 
Since in it t„ is quite generally the surface shearing stress, 
equation (10. 36) must apply in the same way to turbulent flow, too. We 
shall come back to that later. For the special case of vanishing pressure 

drop — sO, equation (10. 36) is transformed into equation (10.4) found 
dx 

before for the flat plate in longitudinal flow. 

The further calculation of the boundary layer on the basis of 
equation (10. 36) is performed for the laminar case according to the method 
of Pohlhausen (reference 24) and for the turbulent case according to the 
method of Gruschwitz (reference 34) (chapter XVIII). 


c. Calculation of the Boundary Layer According to the 


Method of Karman— Pohlhaus en-Holst e in 

For further calculation it is of importance to find a suitable 
expression for the velocity distribution in the boundary layer u(x, y) . 
According to our understanding of the exact solutions of the differential 
equations of the boundary layer this expression must at least satisfy 
the conditions that for y = 0: u = 0, and for y = °°: u = U. Further- 
more the derivative <3u/6y must vanish for large y. Moreover, velocity 
profiles with and without inflection points must be possible, as they 

occur in the pressure decrease and pressure increase region, respectively. 

/ ~ \ 

Finally, a profile with ( ^ ) =0 must be possible in order to have a 

VWy, o 

separation point result from the approximate calculation. 


One chooses for the velocity distribution an expression of the form 
u y ) = Uf(y/5p), and sets, according to Pohlhausen, 
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for f(y/&p) a polynomial of the fourth order, hence: 

4=an + li) 2 + CTi^ + d T) = (10.37) 


valid for 0 — y — 5pi 0 — T) — 1. 6p(x) stands for the boundary layer 

thickness, the dependence of which on x has yet to he calculated. The 
boundary layer thickness of the approximate calculation 5p is here 

provided with the index P ( = Pohlhausen) in order to avoid confusion 
with the boundary layer thickness 5 used before. Whereas for the 
exact solutions the velocity in the boundary layer asymptotically 
approaches the velocity of the potential flow, (u — >U for y — >00 ); 
the value u = U is to be attained in the approximation at a finite 
distance from the surface y = 6p, for reasons of calculation. This 

modification of the actual relation is physically insignificant. 

For the determination of the free constants a, b, c, d in 
equation (10. 37) the following boundary conditions are prescribed, all 
of which follow from the differential equation of the boundary layer 
(equation (8.3)): 


y = 5p: 

y = 0: 


u = U; 


u = 0; 


= 0 . 

d 2 u _ 

8y 

1 

OJ 

& 

1 8 ^u 

u® 

i 

OJ 

& 

dx 


■\ 


> 




(10.38) 


Since the condition of no slip u = 0 for y = 0 is automatically 
satisfied by expression (10.37)* the four free constants a, b, c, d 
are sufficient to satisfy the remaining four conditions. The last 
of the five conditions follows immediately from the exact differential 
equation of the boundary layer if one puts y = 0 and takes the 
boundary conditions into consideration. This condition is particularly 
important since it determines the curvature of the velocity profiles 
near the surface and assures that boundary layer profiles do not acquire 
an inflection point in the region of pressure decrease and do acquire 
one in the region of pressure increase, as required by the exact solution 
according to chapter VXII. From equation (10. 38) follows for the 
coefficients a, b, c, d the equation system: 
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8, + b + c + d = 1 

a + 2b + 3c + 4dL = 0 

2b + 6c + 124 = 0 




> 


(10.39) 


2vUb 


= - U U' 


From the last equation follows immediately: 

2 


b = - i 


5 , 


U* = - 


k 

2 


(10.40) 


4f one introduces the simplification: 


V dx 


(10.41) 


The dimensionless quantity X plays the role of a form parameter of the 
■velocity profiles as will become clear presently. For the remaining 
coefficients one obtains from equation (10.39): 


a = 





(10.40b) 


Hence the expression for the velocity distribution, which satisfies all 
boundary conditions according to equation (10. 37)* reads: 

“ = (2f) - 2 t) 3 + t/*) + | (t) - 3t) 2 + 3t) 3 - t/ 4 ”) 

(10.42) 

g = r(ti) + X g(ti) 

in which 

F(n) = 2 tj - 2t]3 + t) 4 
G (n) = £ (tj - 3t) 2 + 3t] 3 - t/*) 


> 


(10.43) 
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Due to t) = y/Sp(x) the houndary layer thickness 6 p(x) is here the 
only unknown. If that is calculated, the parameter X. follows immediately 
from equation ( 10.41). From equation (10. 42) one understands that tlje 
velocity profiles form with the form parameter A.(x) a one— parameter 
family. The functions F(q) and G(q) indicated in figure 45 and 
table 4 have a universal character, that is, they do not depend on the 
special body shape. The velocity profiles for various values of X are 
plotted in figure 46. The profile with X = 0 is obtained for du/dx = 0, 
that is, for the boundary layer without pressure gradient (flat plate in 
longitidinal flow). The separation profile with (du/cty^ = 0, that is, 

a = 0, has according to equation (10.40b) the parameter X = —12. The 
profile at the stagnation point has, as will be shown below, X = 7.052. 

For X > 12 there result values of u/U >1 in the boundary layer, which 
physically does not make sense. These values therefore have to be excluded. 
Since behind the separation point the boundary layer calculation loses its 
validity anyway, the form parameter X is limited to the region 

- 12 - X - + 12 (10.44) 


The unknown boundary layer thickness 5p(x) remains to be 

calculated. For this the momentum equation (10. 36 ), so far not utilized, is 
at our disposal. Before performing this computation, a few preparatory 
calculations are required, namely the determination of the boundary layer 
characteristics, displacement thickness 5*, momentum loss thickness d, 
and surface shearing stress t q on the basis of the approximation- 

expression equation (10.42): One obtains from equations (10. 32) and 

(10.42): 


5* 


[l ~ F(q) - 4G( t) )J dq 


' q=0 


d_ 

S_ 


[F(q) + 4G(q)] [l - F(q) - XG(q)J dq 


Jq=0 


The calculation of the definite integrals, with the values of F(q) and 
G(q) according to equation (10.43), gives: 


8 * = J X_ 

8 p 10 120 


(10.45) 
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i_ = _37 \ X . 2 

5p 315 9^5 9072 


(10.46) 


Further, there results for the surface shearing stress from t q 
according to equations (10.4-2) and (10.43) 

^0 S _ 12 + X 

» U ' 6 



(10.47) 


Now the momentum equation 
Sp(x). After multiplication by 


(IO. 36 ) is to he used for calculation of 
d/vU it acquires the dimensionless form: 


th3d» 6 *\ U»d 2 _ T o fl 

V \ -8 / V '7u 


(10.48) 


The boundary layer thickness 6 p does not even appear in this equation; 

however, this is not particularly astonishing, since 5p is a rather 

arbitrary quantity of our approximate calculation and therefore without 
special physical significance. The physically important quantities, 
displacement thickness 5* and momentum loss thickness d, appear 
instead in equation (10.48). Hence it suggests itself to first 
calculate from the momentum equation (10.48) and then to pass on 
to 5p by means of equation (10.46). For this purpose one introduces 

according to Holstein and Bohlen (reference 25) aside from the first form 
parameter \ according to equation (10.41) a second form parameter K, 
formed analogously with the momentum thickness 


K 


2 

= U' 
V 


Then one sets: 



V 


(10.49) 


(10.50) 


Then 


K = Z U* 


(10.51) 
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Between the second form parameter K and the first form par ame ter X 
exists, according to equations (10.40) and (10.41) the universal relation: 


-2i _ -x x 

.315 945 9072/ 


(10.52) 


Furthermore, for simplification, one substitutes: 

= 12 — i§o — -=f 1 ( K ) (10.53) 

^7 1 

315 945 9072 


To ^ = To 6p = / 12 + 37 _ _X_ 

p U n U 5p \ 6 /\315 945 



(10.54) 


By introduction of k and Z according to equations (10.49) and. (10. 50) 
and by substituting from equation (10.53) and (10.54) one now obtains 

from the momentum equation (10.48), because of dd*/v - ^ ^ 


| u g + [2 + f i(ic)jK - f 2 (K) = 0 (10.55) 

Finally, one sets as further simplification 

~ 4k - 2^00 K = F(k) (10.56) 


or written in detail: 


F(K) = 2-|l- 


315 9^5 9072 , 


2 - 


116 

315 


X + 


( 2 

V945 120. 


i.) 

20 / 


2 

X + 


9072 


x W10. 57) 


and thus obtains the following differential equation for the momentum 
thickness 


dZ _ F(jO k = 7i U* 
dx TJ 


(10.58) 
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This is a non— linear differential equation of the first order 
for z = $2/^' The fact that the function F(k) is rather complicated 
does not constitute any appreciable drawback, since F(k) is universal, 
that is, independent of the shape of the bc>dy, and hence may be tabulated 
once and for all. The functions F(k), f-jJ*), f 2 ( K ) as also k = k(x) 
according to equation (10.52) are given in table 5* 

As to the solution of equation (10. 58) the following remains to be 
said: The calculation has to start at the stagnation point x = 0. There 

U = 0; and the initial slope ^ would be infinite if F(k) were not 

also equal to zero at the stagnation point. The function F(k) actually 
has a zero which yields a physically significant initial value. This 
zero of F(k) is given when the second bracket in equation (10.57) 
disappears. One finds: 


F(k) =0 for K = k q = 0.0770; X = X = 7.052 (10.59) 


The value X = 7*052 therefore gives the value of the first form 
parameter at the stagnation point. Then the initial slope of the integral 
curve at the stagnation point now has the indeterminate value O, The 

latter may, however, be calculated and hence finally yields the initial 
value and the initial slope of the integral curve as: 



(10.60) 


The index o denotes the values at the stagnation point. With these 
initial values one succeeds easily in performing the integration of 
equation (10. 58), for instance, according to the isocline method. A 
calculation example is given in the appendix, figure 4-7, and table 6. 

The calculation is to be carried up to the separation point X = —12: 
ic = — O.I567. Quantities entering the calculation that are given by the 
potential flow are the velocity U(x) and its first derivative with 
respect to the arc length dU/dx.* (Only at the stagnation point is 
d^rj/dx^ also required for the initial slope of the integration curve.) 

* In PohlhauSen’s treatment (reference 24 ) a differential, equation is 
obtained instead of equation (10. 58), for the quantity z = Sp yv, formed 
analogously to Z. Pohlhausen* s differential equation also contains 
d^U/dx^ which often can be obtained from the given potential flow only 
by a double graphical differentiation. The representation of Holstein 
which completely avoids the quantity d^U/dx^ therefore means an 
essential improvement of the method. 
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The values of the form parameters for the three special cases: 
stagnation point, velocity maximum (U* = 0) (pressure gradient equals 
zero), and separation point are compiled in the table below. 


Case 


K 

Stagnation point 

7.052 

0.0770 

Velocity maximum 

0 

0 

Separation point 

-12 

-0.1567 


The entire process of calculation takes the following course: 

1. The integration of equation (10. 58) yields Z(x), k(x) and 
according to equation (10. 50), also 8(x)j furthermore it yields the 
position- of the separation point. 

2. First form parameter \(x) from equation (10. 52) 

3. Displacement thickness 8* from equation (10.53) 

4. Surface shearing stress t 0 from equation (10.5*0 

5. Boundary layer thickness 8p(x) from equation (10.45) 

6. Velocity distribution u/CJ from equation (10.42) 

Flat Plate in Longitudinal flow 


The special case of the flat plate in longitudinal flow which was 
treated in chapter Xa with a different form for the approximation can 
also be obtained very simply from the present calculation. XJ = U ; 

U t s 0, and hence k = X =0,' and according to equation (10. 58): 


dZ = F(0) = 0.4698 


dx IT 




With the initial value 8=0 for x = 0 there results 


•8 = 0.4698 ^ or 8 = 0.685 W^- 
o v o 


(10.61) 
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whereas the exact value according to Blasius* calculation, equation (9.21h), 
is $ = 0 . 664 \jv x /U Q . Furthermore, it follows that the displacement 
thickness is with b*/$ = f^(0) = 2.54 


S* 



(10.62) 


The shearing stress "becomes, from equation (10.54) with f^O) = 0.235> 


T 


o 


0.343 



(10.63) 


while the exact value t q is t d = 0.332 m.U q \j U 0 /vx according to 
equation (9.17). The agreement with the exact values is rather satisfactory. 
In figure 50b the velocity distribution obtained by the approximate calcula- 
tion also is compared with the exact calculation in the plot u/CJ against 
y/5*. This agreement also is rather good. 


Stagnation Point Profile 


A similar comparison can be performed for the stagnation point 
profile the exact solution of which was given in chapter VI. For this 
case X = A. 0 = 7*052, k = k q = 0.0770. For the approximate calculation 
the momentum thickness is: 


'fo - °- 278 


(10.64) 


wh ereas according to the exact calculation equation (6.34) it is 
$ \JrT t /v =0.292. The displacement thickness is, for the approximate 
calculation; 


5*\H^= fi( K o) \/S = °- 61 ^ 


(10.65) 


vhereas according to th.6 


\ W . 


00^ 


it is 


5* = 0.648. Finally , according to the approximate calculation the 

shearing stress is: 
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To. \IZ = f 2 K o = 0.332 = llq 

nu-\|u« 0.278 


( 10 . 66 ) 


compared with — \ — 

* pU V U* 

the agreement with the 


= cp"(0) = 1.234 by the exact calculation, 
exact values is here also satisfactory. 


Thus 


The velocity distribution of the approximate calculation is compared 
with the exact calculation in figure 50— a in the plot u/U against y/8*. 
Here also the agreement is good. 


The following table contains a compilation of the comparisons Just 
given, between the characteristics 5*, ■Q , t q from exact and from 
approximate calculation. 



Blasius Profile 

Stagnation point profile 

6* \F° 

Vvx 

TO 

5* 

5* JvL 


T o 

6 X 

Vvx 


6 

5 ylr 

1 

d 

n \/u ( 0 3 /v* 

4 u\/v/u» 

Pohlhausen 

Approximation 

1.75 

0.685 

0.343 

2.55 

0.641 

0.278 

1.19 

1 

2.31 

Exact solution 

1.73 

0.664 

0.332 

2.61! 

1 

0.648 

CVJ 

d 

1.234 

2.21 


Of course, it can not yet be concluded from this good agreement of 
the approximate with the exact solution that similar good agreement would 
exist for all the boundary layer profiles along the body. Accurate 
comparisons are not easily performed since very few exact solutions 
reaching from the stagnation point to the separation point exist. However, 
one may conclude from occasionally made comparisons that in the region of 
pressure decrease the agreement is mostly rather good; in the region of 
pressure increase, particularly near the separation point. Borne deviation 
might occur. 

Since no other serviceable methods for boundary layer calculation 
have so far become known, the Pohlhausen method is for the present to be 
regarded as the best. The time required for a boundary layer calculation 
for one side of the body immersed in a given potential flow amounts to 
about three hours ; 

The calculation described here for the plane flow was applied by 
Tomotika (reference 26) to the rotati onally— symmetrical case. 
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d. Examples 

A few examples of "boundary layer calculations will be given, all 
of which were performed according to the approximation method described 
in the previous section. 

The first example, taken from a treatise of H. Schlichting and 
A. Ulrich (reference 35) concerns a series of elliptical cylinders in a 
flow parallel to the major axis. The axis ratios are a^/b^ = 1 

(circular cylinder), 2, 4, 8, 00 (flat plate). The velocity distributions 
of the potential flow are given in figure 51* For slender cylinders they 
are characterized by the fact that the velocity has a very flat maximum 
and decreases steeply toward the front and rear. The result of the 
boundary layer calculation is illustrated in figure 52. Figure 52-a 
shows the dimensionless displacement thickness as a function of the arc 
length. It is noteworthy that for all elliptic cylinders the boundary 
layer on the front half is smaller than for the flat plate; this is 
caused by the pressure decrease. The laminar— separation point for the 
circular cylinder lies at s/t = 0.604 (9 = 109.5°)» With increasing 

fineness ratio it shifts further toward the rear; it is also plotted in 
figure 51* Figure 52-b shows the variation of the form parameter X, and 
figure 52— c the variation of the surface shearing stress. For every 
elliptic cylinder the latter has a maximum, the position of which shifts 
frontward with increasing fineness ratio. The variation of the shearing 
stress for the cylinder of axis ratio eight shows only an insignificant 
difference from the one for the flat plate in longitudinal flow. Figure 53 
gives a survey of the variation of the boundary layer on the body and the 
velocity profiles at various points along the surface. Corresponding 
calculations for the rotationally— symmetrical case (that is, for ellipsoids 
of revolution with flow parallel to the axis) have been performed by 
Pretsch (reference 27). 

The second example (reference 35) gives the boundary layer on a 
symmetrical Joukowski profile of fifteen percent thickness for lift 
coefficients in the region of c a = 0 to 1. Figure 5^ gives the velocity 

distribution according to ,the potential theory for the lift coefficients 
c a = 0, 0.25, 0.50, 0.75, and 1.0. For the symmetrical approach flow 

(c a = 0) the velocity maximum lies rather far toward the front at 

s/t = 0.141. With growing lift coefficient the velocity maximum increases 
on the suction side and decreases on the pressure side. Simultaneously 
the maximum shifts farther forward on the suction side and farther rearward 
on the pressure side. The magnitude of the velocity maximum and its 
position are of primary importance for the development of the boundary 
layer and in particular for the position of the separation point. The 
results of the boundary layer calculation are plotted separately for the 

suction and pressure sides in figure 55. Figure 55-a gives the 
variation of the displacement thickness, figure 55-b the form • 

parameter \,,and figure 55-c the surface shearing stress. The position 
of the laminar separation point is plotted in figure 5^-. Figure 5 6 gives 
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a survey of the variation of the "boundary layer along the surface and of 
the velocity distribution in the boundary layer for c a = 0, figure 57 
for c a = 1. 

Finally, the third example gives a survey- of the influence of the 
most important profile parameters of a wing profile on the laminar 
boundary layer. K. Bussmann (reference 36) performed the boundary layer 
calculation for a family of Joukowsky profiles of relative thickness 
d/t = 0 to 0.25 and relative camber f/t = 0 to 0.08 for c a — values 

from 0 to 1. Figure 58 shows the family of Joukowsky profiles. Of the 
very voluminous results, only the position of the separation point shall 
be shown here. Figure 59 shows the position of the separation point on 
the suction side as a function of thickness, camber, and lift coefficient; 
figure 60 shows the same result for the pressure side. 

Herewith the discussion of the approximation method for calculation 
of the laminar boundary layer will be concluded. 


Eighth Lecture (January 2 6, 1942) 
CHAPTER XI. PREVENTION OF SEPARATION 


For practical flow problems the flow with pressure increase 
(retarded flow) plays an important role. It is always desirable that 
no separation o f the flow from the wall occur, because of the resulting 
large losses in energy. The wing presents a good example. A pressure 
increase exists on the suction side toward the trailing edge (fig. 6l). 

If separation occurs, the wing will have an undesirably large drag and 
small lift. Another example is the flow in an expanding passage 
(diffuser) which transforms kinetic energy into pressure energy (as for 
instance in the wind tunnel or in the bucket grid of a turbine). 

Calculations will presently show that the ability of a laminar flow 
to overcome a pressure increase without separation is exceedingly small. 
Thus the pressure increases present in practical flows would, for laminar 
flow, almost always lead to separation. The reason that, nevertheless, 
in many cases of practical flows considerable pressure increases are 
surmounted without separation is that the flow is turbulent . As we shall 
see more clearly later, the ability to overcome a pressure increase with- 
out separation is very much greater for turbulent than for laminar flow. 
Since, moreover, the pressure increase always gives rise to an early 
transition from laminar to turbulent*, one has to deal almost exclusively 
with turbulent flow in practical flows with pressure increase. 

Nevertheless it is useful to clarify the fundamental relations 
regarding prevention of separation for laminar flow, particularly because 
its phenomena lend themselves more readily to numerical treatment than 
those of the turbulent flow. 


*Details are given in chapter XXI "Transition from Laminar to Turbulent." 
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Various possibilities exist for prevention of separation. The 
simplest way is to make the pressure increase so small that separation 
is avoided. A numerical estimation in the next section will give 
information about this possibility. Another possibility consists in 
artificially influencing the boundary layer, for instance by blowing 
or suction of fluid, or else by application of an auxiliary wing that 
provides acceleration at the critical points of the boundary layer. 
Some details will be given in the following sections of this chapter. 


a. Estimation of the Admissible Pressure Gradient 

We are going to make, following Prandtl (reference 2), a generally 
valid estimation of the pressure increases in a laminar boundary layer 
that are possible without the occurrence of separation. We take as the 
basis the Karmen— Pohlhausen approximate calculation according to 
Chapter X and make the assumption that under the effect of the pressure 
gradient given by the potential flow the boundary layer has developed 
till near the separation point (Point 0 in fig. 62). From here on 
the pressure distribution is to be such that the form of the velocity 
profile does not change further downstream, that is, the form parameter 
X is to remain constant. Since the value corresponding to the separation 
point is X = —12, this constant X— value shall be chosen at X = —10*. 

A definite value of the second form parameter (according to table 5) 
corresponds to this choice: 


X = -10; k = -O.I369; F(k) = 1.523 (11.1) 


For the prevention of separation the following relation between the 
potential— flow velocity U(x) and the momentum thickness d(x) results 
according to equations (10. 50) and (10. 51): 


or 



Z = 


O.I369 

-U* 


( 11 . 2 ) 


dZ 

dbc 


0.1369 



(11.3) 


or 


V H = 0.1369 = 0.1369 <T 


(11.4) 


if one puts for simplification 

* 

• nir " J 

*At any rate, this X-value must be negative, since otherwise the flow 
in question is not retarded. 
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(IX. 5) 


On the other hand, the momentum equation according to 
equation (10. 58) holds for the further development of the boundary layer 
for x > 0: 


U dx = F ( K ) = F (“°" 1 3 6 9) - 1.523 


( 11 . 6 ) 


The numerical value must be substituted for F(k ) , if the form parameters 
are to remain constant at the values given by equation (11.1). From 
equations (11.6) and (11.4) follows therewith, for the constancy of the 
form parameter X = —10, the conditional equation 

O.I369 ~ = 1.523 

U' 


0=22^=11.13^11 (11. 7 ) 

u » 2 


For o > 11 the boundary layer can still bear the pressure increase; 
for a < 11 separation occurs; for 0 = 11 the boundary layer always 
remains with X = —10, on the verge of separation. Qualitatively, the 
following can be immediately said about the distribution of the potential- 
flow velocity U(x) which gives no separation. Because of equation (11. 7) 
a necessary condition for avoiding separation in retarded flow is: 


or 




U” > 0 f or U» < 0 


(11.7a) 


that is, a negative velocity gradient U* must exist, the magnitude of 
which decreases in the flow direction. If, therefore, the curve U(x) 
in figure 63 is curved downward behind the maximum (U* < 0), separation 
occurs in ever y case; if it is curved upward (U ,, > 0), separation some- 
times does not occur. The limiting case U w = 0 for U f < 0 always 
leads to separation. The sufficient condition for avoiding separation 
is UU"AJ‘ 2 > 11. 
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One now proceeds to calculate what potential flow and what boundary 
layer thickness variation correspond to a = +11. From, equation (11.7) 
follows : 


= 11 21 

u> u 


and after integration: log U* = 11 log U — log C£ or 


U» 

ttII 


C* 


with as integration constant. Repeated integration gives: 


1 TT -10 

— U = C* I + Co 

10 1 2 


( 11 . 8 ) 


For x = 0, U(x) shall he U(x) = U 0 , thus 


C — _i_ tt ~1° 

C 2 - 10 U o 


(11-9) 


Furthermore, one puts 


Uo 


10 


= C, 


( 11 . 10 ) 


and obtains from equation (11.8) for the potential flow 


U = 


U„ 


(1 + 10 C.,*:) 


0.1 


( 11 . 11 ) 


Thereby is found the desired velocity distribution that Just avoids 
separation. The constant can be determined from the boundary layer 

thickness S Q at the initial point x = 0: 


lo8 
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According to equation (11.11) 


U* 


C 1 U o 


(1 + 10 C-,1) 


1.1 


and thence 


5 = 


(1 + 10 Cl x )°- 55 

1 o 


From 5 = &p for x = 0 follows 


n - 10 v 

1 “ 

V>o 


( 11 . 12 ) 


and thus, as the final solution for the potential flow and the boundary 
layer thickness variation. 


u = U G (1 + 100 -2S-\ 

U o 5 o ) 


- 0.1 


(11.13) 


6 = 8 Q [1 + 100 — 

W 


0.55 


( 11 . 1*0 


The permissible retardation (velocity decrease) is therefore comparable 
10_ 

to ±/\fx and is thus very small. The velocity is thus very close to 
the constant velocity of the flat plate in longitudinal flow. For the 
present case the growth of the boundary layer thickness 5 must therefore 

1/2 

be somewhat larger than for the flat plate, where 8 ~ x . Here 
0.55 

8 ** x ; thus the increase is only slightly larger. 

The flow in a divergent channel with plane walls (two-dimensional 
problem) will be treated as another example. In figure 6k let x be 
the radial distance from the origin 0. The walls start at x = a, where 
the entrance velocity of the potential flow equals U Q . The potential 
flow is 

♦Compare Chapter IX e where it was found, as exact solution of the 
differential equation of the boundary layer, that in retarded flow 

-0.091 


separation occurs when U(x) = u^ x 
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"(*> = U 0 j 


U» = - U 0 — 

0 2 

■XT 


> 


(11.15) 


U" = 2XJ — 
o 

X-* 




ThuB U* < 0 and U" > 0 for all x bo that the necessary condition 
equation (11.7a) for avoiding separation is satisfied. However, calcula- 
tion of the dimensionless number a according to equation (11.5) gives 


o = 2 


( 11 . 16 ) 


The sufficient condition for avoiding separation, a > 11 according to 
equation (11.7) is therefore violated. For the divergent channel with 
plane walls separation therefore occurs for any included angle. This 
example shows especially clearly the low ability of the laminar flow to 
overcome a pressure increase without separation. According to a calcula- 
tion of Pohlhausen (reference 24) the separation point lies at (x/a) = 1.21 

and thus is independent of the included angle a. 


b. Various Technical Arrangements for Avoiding Separation 

It is a favorable circumstance for technical applications that for 
higher Reynolds numbers the boundary layer does not remain laminar but 
becomes turbulent. The turbulence consists of an irregular mixing motion. 
By this mixing motion momentum is continuously transported into the layers 
near the wall, and the particles retarded at the wall are carried out into 
the free stream and thus re-accelerated. 

Because of this mechanism the turbulent flow is able to withstand, 
without separation, considerably higher pressure increases than the 
laminar flow; thus the pressure increases existing in technical flows 
are made possible. 

A few technical possibilities for avoiding separation will be 
discussed. 

1. Blowing. For a wing profile the separation of the boundary 
layer for large angles of attack (fig. 65 ) can be prevented by blowing 
air in the flow direction from a slot directed toward the rear. The 
velocity for the layer near the surface is thus increased by the energy 
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supplied and the danger of separation is therefore eliminated. It is 
true that in the practical execution not much is gained, because of the 
large jet energy required for any considerable improvement of the flow. 

In order to make the energy output small, the width of the Jet must be 
kept small. But then the Jet, soon after its exit, breaks up into 
vortices. 

2. Another possibility of avoiding separation is the arrangement 
of a slotted wing according to figure 66 . The effect depends on the 
boundary layer formed on the slot AB being carried away into the free 
stream, before it separates, by the flow through the slot. A new 
boundary layer develops at C which is, however, at first very thin 
and reaches D without separation. 

The same principle is used for the Townend ring and NACA cowling 

(fig. 67). 

3 . Suction. A further possibility for the prevention of separation 
is suction. For the wing, for instance, the retarded boundary— layer 
material is. sucked off into the interior of the wing through one or 
several slots (fig. 68 ). The point of suction lies slightly ahead of 

or behind the separation point so that no reversal of the flow can occur. 

A new boundary layer which at first is very thin develops behind the 
suction point and permits the pressure to increase further. In this 
manner one can overcome considerably larger pressure increases and 
attain higher values of maximum lift for the wing. Many different suction 
arrangements for increasing maximum lift have been investigated by 
0. Schrenk (reference 28). Values for c a of 3 to k were obtained. 


c. Theory of the Boundary Layer with Suction 

Suction is a very effective means for influencing the friction 
layer on a body immersed in a flow and particularly for avoiding 
separation. This was pointed out for the first time in 190 k by L. Prandtl 
in his fundamental work on the boundary layer. 

Another possibility of application of suction, recognized only 
recently, is to keep the friction layer laminar. Here the boundary layer 
is, by suction, kept so thin that transition to the turbulent state of 
flow is avoided. The surface friction drag is thereby reduced. Experimental 
investigations of this effect were carried out by Ackeret (reference 39 ). 

The laminar friction layer with suction can also be subjected to a 
numerical treatment which will be briefly discussed. The following 
assumptions are made- for th 9 calculation: 
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1. The suction is introduced into the calculation through the 

assumption that the normal velocity at the wall v 0 (x) is 
different from zero. The wall is therefore assumed to he 
permeable. A continuous distribution of the suction velocity 
v 0 (x) serves the purpose of numerical treatment best. 

2. The suction quantities are so small that only the parts in the 

immediate neighborhood of the wall are sucked from the 
boundary layer. This leads to a very am«.i 1 ratio of suction 
velocity v Q (x) to free stream velocity U 0 : v 0 /U Q = 0.001 

to 0.01. 

3. The no— slip condition at the wall u = 0 is retained with 

suction, likewise the expression for the wall shearing stress 



The equations of motion for the boundary layer with Buction therefore 
read 


8u . du 
u — + v — 

<bc dy 




du | dv 
dx dy 


0 


J 


with the boundary conditions 


(11.17) 


y=0 u=0 v= v Q (x) 

y = oo u = U 

v 0 < 0 signifies suction; v Q > 0 blowing. 





(11.18) 


As in chapter X b the momentum theorem is again applied to the 
boundary layer with suction. The momentum equation for the boundary layer 
with suction is obtained in exactly the same manner as in chapter X b 
(compare fig. MO provided one takes into consideration, in addition, 
that the nor mal velocity at the wall is different from zero. In chapter 
X b the momentum equation was derived by integration of the equation of 
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motion for the x— direction over y between the limits y = 0 and 
y = oo (compare equation (10.28)). One imagines exactly the same calcula- 
tion performed for the boundary layer with suction: then the expression 

for the normal velocity at the distance from the wall y = h is 
different, compared with the calculation in chapter X b. The normal 
velocity now becomes 


v, = 


v — 
o 


R h 


8u Arr 

5 ^ 


do 


(11.19) 


The remaining calculation is exactly the same as in chapter X b and 
finally yields as the momentum equation for the boundary layer with 
suction. 


1° = U 2 — + ( 2-3 + 6*) U — — v n U 
p dx dx 0 


( 11 . 20 ) 


The newly added term — v Q U (compared with equation (10. 36)) gives the 
loss of momentum due to the suction at the wall. 


We shall now treat the special case of the flat plate with suction 
in longitudinal flow (fig. 69) (reference 29). The free-stream velocity 
is U 0 . Equation (11.20) then becomes 


2 






( 11 . 21 ) 


if one takes the law for the laminar wall shearing stress into consider- 
ation. Furthermore, the assumption is made that the suction velocity 
(or blowing velocity) — v 0 along the plate is constant. In this case 
one can obtain from the momentum equation (11.21), by the following 
, simple calculation, an estimate of the variation of the momentum thick- 
ness along the plate. One puts 



( 11 . 22 ) 


p — 0 signifying a dimensionless form parameter of the velocity profile . 
It may be assumed, to a first approximation, that p varies only little 
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with the length of run x; accordingly, 0 will be considered constant. 
Then’ equation (11.21) may be written 


1 dd_ 

2 dx 



(11.23) 


with the initial condition d = 0 for x = 0. For suction (v Q < 0) one 
obtains dd /dx = 0 for 


= P (suction) 


(11.24) 


(that is, therefore, the momentum thickness reaches, after a certain 
approach length, a constant asymptotic value given by equation (11.24)). 
Simultaneously, displacement thickness, velocity distribution, and all 
other boundary layer coefficients also become asymptotically independent 
of x. 


For blowing (v Q > 0) the value dd/dr is, according to equation 

( 11 . 23 ), larger than zero along the entire plate; that is, d(x) increases 
with the length of run x without limit so that for large values of x, 
one can neglect in equation ( 11 . 23 ) the first term on the right side as 
compared with the second. One obtains therefore, as asymptotic law, 

d^ * ^ x (blowing) ( 11 . 25 ) 


On the whole, one obtains the remarkable result that for the flat plate 
in longitudinal flow with constant suction or blowing velocity, the 
boundary layer thickness for suction becomes constant after a certain 
approach length, whereas for blowing, it increases proportionally to the 
length of run x. In between lies the case of the impermeable wadi 
where the boundary layer thickness increases with Vx. 

For the case of the laminar boundary layer with the asymptotically 
constant boundary layer thickness it is also possible to give immediately 
an exact solution of the differential equations of the boundary layer in 
a surprisingly simple form. In this case du/cbc = 0, hence also, 
according to equation (11.17), &v/$y = 0 and therefore 


v(x,y) = v^ = constant 


( 11 . 26 ) 
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Hence there follows from equation (11. 17) 


v 


o 


du 

Sy 



(11.27) 


and from it the solution which satisfies the houndary conditions 
equation ( 11 . 18 ) 


u(y) = u o < l 



( 11 . 28 ) 


v J 

From this equation results the displacement thickness of the asymptotic 
boundary layer 


8 * 


00 



(11.29) 


the momentum thickness 


d 

00 




(11.30) 


and the form parameter 



By comparing equation (11.29) with 


equation (11.24) one finds the factor 0 = 1. The velocity distribution 
of the asymptotic boundary layer profile according to equation ( 11 . 28 ) 
is plotted in figure 70 together with the Blasius solution for the 
impermeable wall. 


Herewith the considerations of boundary layer with suction will 
be concluded. 


CHAPTER XI I. APPENDIX TO PART I 
a. Examples of the Boundary Layer Calculation 
According to the Pohlhausen— Holstein Method 


For the integration of the differential equation (10. 58 ) it is 
best to use the isocline method. It is expedient to calculate with 
dimensionless quantities. The arc length s is made dimensionless by 
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dividing "by a characteristic length of the "body immersed in the flow, 
for instance, for the wing, "by the wing chord t. The variable Z = 

U 0 

is made dimensionless by multiplying by — . Thus one puts: 


Z* 



V t 




x* 


8 


= 7 


J 


( 12 . 1 ) 


Hence the differential equation reads: 


— — — — F( ic) j R =±™ Z * 
dx* U ds 


( 12 . 2 ) 


The calculated example concerns a symmetrical wing profile (j 015) in 
symmetrical approach flow (c a = 0). The prescribed potential— flow 
velocity and its first derivative with respect to the arc length is 
given in tpble 6. The initial values for the integration are calculated, 
according to equation (10.60), to be, for the present case: 


Z 0 * =0.00 149 



t 


since at the stagnation point d^U/ds^ = 0. The auxiliary function F( k) 
required for the integration is given in figure 47-a and table 5* The 
calculation according to the isocline method is shown in figure 48. Here 
the curve k = -O .1567 which gives the separation point can be calculated 
according to the relation: 


z * = _k/_L £2 = _ 0.1567/ — 

A/ U 0 ds ‘ " 


t dU 


U 0 ds 


The intersection of the integral curve with this curve gives the separa- 
tion point. As a result of the integration one obtains at first the varia- 
tion of the momentum thickness. By means of the function b*/$ = f( K ) 


and 


3 . 1 ° 
u n 


= f 2 (i0 


given in table 5 one can also calculate the displacement 


n6 
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thickness and the shearing stress. The result of the calculation is 
compiled in table 6 and given in figure h9- Moreover, the velocity 
distribution in the boundary layer can be seen from figure 56 . 


Translated by Mary L. Mahler 
National Advisory Committee 
for Aeronautics 
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TABLE I. - THE FUNCTION cp OF THE PLANE STAGNATION POINT 
FLOW (ACCORDING TO HIEMENZ (REFERENCE 10)); TO FIGURE l6 


1 

9 

cp « 

cp ” 

0 

0 

0 

1.23264 

0.1 

0.0060 

0.1183 

1.1328 

0.2 

0.0233 

0.2266 

1.0345 

0.3 

0.0510 

0.3252 

0.9386 

0.4 

0.0881 

0.4144 

0.8463 

0.5 

0.1336 

0.4946 

0.7583 

0.6 

0.1867 

0.5662 

0.6751 

0.7 

0.2466 

0.6298 

0.5973 

0.8 

0.3124 

0.6859 

0.5251 

0.9 

0.3835 

0.7350 

0.4586 

1.0 

0.4592 

0.7778 

0.3980 

1.1 

0.5389 

0.8149 

0.3431 

1.2 

0.6220 

0.8467 

0.2937 

1.3 

O.7081 

0.8739 

0.2498 

1.4 

0.7966 

0.8968 

0.2109 

1-5 

0.8873 

0.9161 

0.1769 

1.6 

0.9798 

0.9324 

0.1473 

1.7 

1.0738 

0.9457 

0.1218 

1.8 

1.1688 

0.9569 

0.0999 

1.9 

1.2650 

0.9659 

0.0814 

2.0 

1.3619 

0.9732 

0.0658 

2.1 

1.4596 

0.9792 

0.0528 

2.2 

1-5577 

0.9841 

0.0420 

2.3 

1.6563 

0.9876 

0.0332 

2.4 

1.7552 

O.9905 

0.0260 

2.5 

1.8543 

0.9928 

0.0202 

2.6 

1.9537 

0.9946 

0.0156 

2.7 

2. P 533 

0.9960 

0.0119 

2.8 

2.1529 

0.9971 

0.0091 

2.9 

2.2528 

0.9979 

0.0068 

3.0 

2.3525 

0.9985 

0.0051 

3-1 

2.4523 

0.9988 

O.OO36 

3-2 

2.5522 

0.9992 

0.0027 

3.3 

2.6521 

0.9994 

0.0023 

3-4 

2,7521 

0.9996 

0.0019 

3-5 

2.8520 

0.9997 

0.0014 

3.6 

2.9520 

O.9998 

0.0010 

3-7 

3.0519 

0.9999 

0.0008 

3.8 

3.1518 

0.9999 

0.0004 

3-9 

3.2518 

0.9999 

0.0003 

4.0 

3-3518 

1.0000 

0.0002 

4.1 

3.4518 

1,0000 

0.0001 

4.2 

3.5518 

1.0000 

0.0001 

4.3 

3.6518 

1.0000 

0 
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TABLE H. - THE FUNCTION f OF THE BOUNDARY LAYER ON THE FLAT ELATE 
IN LONGITUDINAL FLOW (ACCORDING TO BLASHJS (REFERENCE 8)); TO 


FIGURE 30 tj = yl| = f*(ij) 

1 VI U 0 


n 

f 

f* 

f" 

0 

0 

0 

0.33206 

0.2 

0.0066I 

0.06611 

0.33199 

0.1 

0.02656 

0.13277 

0.33117 

o .6 

0.0597I 

O.I989I 

0.33008 

0.8 

0.10611 

0.26171 

0.32739 

1.0 

0.16557 

0.32979 

0.32301 

1.2 

0.23795 

0.39378 

0.31659 

i.i 

0.32298 

0.45627 

0.30787 

1.6 

0.I2032 

O.51676 

0.29917 

1.8 

0.52952 

0.57177 

0.28293 

2.0 

0.65003 

0.62977 

0.26675 

2.2 

0.78120 

0.68132 

0.21835 

2.1 

0.92230 

0.72899 

0.22809 

2.6 

1.07252 

0. 77216 

0.20616 

2.8 

1.23099 

0.81152 

0.18101 

3.0 

1.39682 

0.81605 

0.16136 

3.2 

1.56911 

0.87609 

0.13913 

3.4 

1.7I696 

0.90177 

0.11788 

3.6 

1.92954 

O.92333 

0.09809 

3.8 

2.11605 

0.9U12 

0.08013 

1.0 

2.30576 

0.95552 

0.06121 

1.2 

2.19806 

0.96696 

0.05052 

1.1 

2.69238 

0.97587 

0.03897 

1.6 

2.88826 

0.98269 

0.02918 

1.8 

3.08531 

0.98779 

0.02187 

5.0 

3.28329 

0.99155 

0.01591 

5.2 

3.I8189 

0.99I25 

0.01131 

5.1 

3 . 68091 

0.99616 

0.00793 

5.6 

3.88031 

0.99718 

0.00513 

5-8 

I.07990 

0.99838 

0.00365 

6.0 

I.2796I 

0.99898 

0.00210 

6.2 

1.17918 

0.99937 

0.00155 

6.1 

1.67938 

0.99961 

0.00098 

6.6 

1.87931 

0.99977 

0.00061 

6.8 

5.07928 

0.99987 

0.00037 

7.0 

5.27926 

0.99992 

0.00022 

7.2 

5.17925 

0.99996 

0.00013 

7.1 

5.6792I 

0.99998 

0.00007 

7.6 

5.87921 

0.99999 

o.ooool 

7.8 

6.07923 

1.00000 

0.00002 

8.0 

6.27923 

1.00000 

0.00001 

8.2 

6.17923 

1.00000 

0.00001 

8.1 

6.67923 

1.00000 

0 

8.6 

6.87923 

1.00000 

0 

8.8 

1 

7.07923 

1.00000 

0 
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TABLE III. - THE FUNCTIONS f , g AND h OF THE BOUNDARY LAYER 

3 5 5 

ON THE CYLINDER (SYMMETRICAL CASE) ACCORDING TO HO WART H 
(REFERENCE 15) AND FROSSLING (REFERENCE l6). n = y 



f 3 

f 

3 

f 3 " 


g 5 

g 5 

g" 

8 5 

0 

0 

0.7244 

0 

0 

0 

0.6348 

0.0035 

0.0675 

0.6249 

0.2 

0.0114 

0.1072 

0.4402 

0.0132 

0.1251 

0.5286 

O.k 

0.0405 

0.1778 

0.2717 

0.0282 

0.1734 

0.4375 

0.6 

0.0806 

0.2184 

0.1408 

0.0476 

0.2129 

0.3539 

0.8 

0.1264 

0.2367 

0.0483 

0.0705 

0.2444 

0.2780 

1.0 

0.1742 

0.2399 

-0.0106 

0.0962 

0.2688 

0.2112 

1.2 

0.2218 

0.2342 

-0.0431 

0.1240 

0.2869 

0.1530 

l.k 

0.2676 

0.2239 

-O.O567 

0.1534 

0.2997 

0.1037 

1.6 

0.3112 

0.2123 

-O.O58O 

0.1838 

0.3080 

0.0626 

1.8 

O.3526 

0.2012 

-0.0522 

0.2149 

0.3125 

0.0292 

2.0 

O.3918 

0.1916 

-0.0432 

0.2462 

0.3140 

0.0028 

2.2 

0.4293 

0.1839 

-0.0335 

0.2776 

0.3132 

-0.0173 

2 .k 

0.4655 

0.1781 

-0.0245 

0.3088 

0.3107 

-0.0320 

2.6 

0.5007 

0.1740 

-0.0171 

0.3397 

O.3070 

-0.0420 

2.8 

0.5352 

0.1712 

-0.0114 

0.3702 

O.3025 

-0.0482 

3.0 

O.5692 

O.1694 

-0.0072 

0.4002 

0.2947 

-0.0513 

3-2 

0.6030 

0.1682 

-0.0043 

0.4297 

O.2923 

-0.0518 

3-4 

O.6365 

O.I676 

-0.0026 

0.4587 

0.2871 

-0.0506 

3.6 

0.6700 

O.1672 

-0.0015 

0.4871 

0.2822 

-0.0480 

3.8 

0.7034 

O.1669 

-0.0010 

0.5151 

0.2775 

-0.0444 

4.0 

0.7368 

0.1668 

-0.0004 

0.5426 

0.2733 

-0.0402 

4.2 

0.7701 

0.1667 

-0.0001 

0.5698 

0.2695 

-O.O358 

4.4 

0.8035 

0.1667 

-0.0001 

O.5966 

0.2662 

-0.0314 





0.6230 

0.2632 

-O.O27I 


K 

K 

K 

0.6492 

0.2607 

-0.0230 


5 

5 

5 

0.6752 

0.2586 

-0.0194 

0 

0 

0 

0.1192 

0.7010 

0.2568 

-0.0160 

0.2 

0.0017 

0.0141 

0.0249 

0.7266 

0.2554 

-0.0131 

0.4 

0.0045 

0.0117 

-0.0436 

0.7520 

0.2542 

-0.0106 

0.6 

0.0057 

-0.0010 

-0.0783 

0.7774 

0.2533 

-O.OO85 

0.8 

0.0039 ! 

-0.0176 

-0.0833 

0.8027 

O.2525 

-0.0067 

1.0 

-0.0012 

-0.0330 

-O.O68O 

0.8279 

O.2519 

-0.0052 

1.2 

-0.0090 

-0.0441 

-0.0423 

0.8531 

0.2515 

-0.0041 

1.4 

-0.0185 

-0.0498 

-0.0149 

0.8782 

0.2511 

-0.0032 

1.6 

-0.0286 

-0.0503 

+0.0088 

O.9033 

0.2508 

-0.0024 

1.8 

-0.0384 

-0.0468 

O.O256 

0.9284 

0.2506 

-0.0019 

2.0 

-0.0472 

-0.0406 

0.0351 

0.9534 

0.2504 

-0.0014 

2.2 

-0.0546 

-O.O331 

O.O38O 

0.9785 

0.2503 

-0.0011 

2.4 

-0.0604 

-O.0257 

O.O36I 

I.OO35 

0.2502 

-0.0008 

2.6 

-0.0649 

-O.OI89 

0.0312 

1.0285 

0.2502 

-0.0006 

2.8 

-0.0681 

-O.OI33 

0.0249 

1.0535 

O.25OI 

-0.0004 

3.0 

-0.0703 

-O.OO89 

0.0187 

1.0785 

O.25OI 

-0.0003 

3.2 

-0.0717 

-O.OO58 

O.OI32 

I.IO35 

O.25OO 

-0.0002 

3.4 

-0.0726 

-O.OO36 

0.0089 

1.1285 

O.25OO 

-0.0001 

3.6 

-0.0732 

-0.0021 

0.0057 




3.8 

-0.0735 

-0.0012 

O.OO36 




4.0 

-0.0737 

-0.0006 

0.0022 




4.2 

-0.0738 

-0.0003 

0.0012 


1 


4.4 

-0.0738 

-0.0001 

0.0007 
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4 


TABLE 17. - THE FUNCTIONS F(y/Sp) and G(y/6 p ) FOE THE TELOCITY 
DISTRIBUTION IN TEE BOUNDARY LAYER ACCORDING TO POHLHAUSEN 
(REFERENCE 24) AND HOWAETH (REFERENCE 15) 


7 

5P 

F 

G 

0 

0 

0 

0.1 

0.1981 

0.01215 

0.2 

O.3856 

0.01725 

0.3 

0.5541 

0.01715 

0.4 

0.6976 

0.0144 

0.5 

0.8125 

0.0104 

0.6 

0.8976 

0.0064 

0.7 

0.9541 

0.00315 

0.8 

O.9856 

0.00105 

0.9 

O.998I 

0.00015 

1.0 

1 

0 
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TABLE V. - AUXILIARY FUNCTIONS FOR THE BOUNDARY LAYER CALCULATION 
ACCORDING TO HOLSTEIN (REFERENCE 25 ) 


X 

K 

F(k) 

- f 

f 2 (‘) 

* U H 

15 

0.0885 

-0.0657 

2.279 

0.345 

14 

0.0920 

- 0.0814 

2.262 

0.351 

13 

0.0941 

-O.O913 

2.253 

0.354 

12 

0. 0911-8 

-0.0946 

2.250 

0.356 

11 

0.09U1 

- 0.0911 

2.253 

0.354 

10 

0.0920 

- 0.0806 

2.260 

0.351 

9 

0.0882 

-0.0608 

2.273 

0.346 

8 

0.0831 

-0.0332 

2.289 

0.340 

7.8 

0.0820 

- 0.0271 

2.293 

0.338 

7.6 

0.0807 

-0.0203 

2.297 

0.337 

7.4 

0.0794 

- 0.0132 

2.301 

0.335 

7.2 

0.0780 

- 0.0051 

2.305 

0.333 

7.052 

0.0770 

0 

2.308 

0.332 

7 

0.0767 

0.0021 

2.309 

0.331 

6.9 

0.0760 

0.0061 

2.312 

0.330 

6.8 

0.0752 

0.0102 

2.314 

0.330 

6.7 

0.0744 

0.0144 

2.316 

0.329 

6.6 

0.0737 

0.0186 

2.318 

0.328 

6.5 

0.0729 

0.0230 

2.321 

0.327 

6.4 

0.0721 

0.0274 

2.323 

0.326 

6.3 

0.0713 

0.0319 

2.326 

0.325 

6.2 

0.0706 

0.0365 

2.328 

0.324 

6.1 

0.0697 

C .0412 

2.331 

0.322 

6 

0.0689 

0 - 0459 

2.333 

0.321 

5 

0.0599 

0.0978 

2.361 

0.310 

4 

0.0497 

0.1579 

2.392 

O.297 

3 

0.0385 

0.2255 

2.427 

0.283 

2 

0.0264 

0.3000 

2.466 

0.268 

1 

0.0135 

0.3820 

2.508 

0.252 

0 

0 

0.4698 

2.554 

O.235 

-1 

- 0.0140 

0.5633 

2.604 

0 .?L 7 

—2 

- 0.0284 

0.6616 

2.658 

0.199 

-3 

-0.0429 

0.7640 

2.716 

0.179 

-4 

-0.0575 

0.8698 

2.779 

0.160 

-5 

-0.0720 

0.9780 

2.847 

0.140 

-6 

-0.0862 

1.0853 

2.921 

0.119 

-7 

-0.0999 

1.1981 

2.999 

0.100 

-8 

-0.1130 

1.3078 

3.084 

0.079 

-9 

-0.1255 

1.4173 

3.177 

0.059 

-10 

-0.1369 

1.5231 

3.276 

0.039 

-11 

-0.1474 

1.6251 

3.383 

0.019 

-12 

-0.1567 

1.7237 

3.500 

0 

-13 

- 0.1648 

1.8159 

3.627 

-0.019 

-14 

- 0.1715 

1.9020 

3.765 

-0.037 

-15 

- 0.1767 

1.9821 

3.920 

-0.054 






TABLE VI. - EXAMPLE FOR THE BOUNDARY LAYER CALCULATION ACCORDING 
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Figure 4.- The shearing stress (to fig. 3). 


y 



Figure 5.- The deformation of a pure elongation. 



Figure 6.- Pure angular deformation (e > 0). 



128 


NACA TM No. 1217 


Figure 8.- 


Figure 9.- 


y 



Figure 7.- Pure angular deformation (f > 0). 



Analogy between heat boundary layer and flow boundary layer 


X 



Types of solutions of the Navier -Stokes differential equations 
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Figure 16.- The velocity profile of the plane stagnation point flow. 



Figure 17.- The convergent and divergent channel. 



Figure 18.- Velocity distribution in the convergent channel. 
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77777777 ' 
Reverse flow 


Figure 19.- Velocity distribution in the divergent channel. 



Figure 21.- Viscous flow around a sphere. 
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Figure 22.- Streamline pattern of the viscous flow around a sphere 

(according to Stokes). 



Figure 23.- Streamline pattern of the viscous flow around a sphere 

(according to Oseen). 


u 



Figure 24.- Concerning Prandtl’s boundary-layer equation. (Boundary - 

layer thickness S magnified.) 
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Figure 28.- 


Concerning the calculation of the friction drag. 
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Figure 29.- The boundary layer on the flat plate in longitudinal flow. 



Figure 30.- Velocity distribution u(x,y) in the boundary layer on the 

flat plate (according to Blasius). 



Figure 31.- The transverse velocity v(x,y) in the boundary layer on 

the flat plate. 


Figure 32.- 



The boundary layer on a cylindrical body of arbitrary 
cross section (symmetrical case). 
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Figure 33.- The function fg' of the velocity distribution in the 

boundary layer. 



Figure 34.- The functions g§' and h^' of the velocity distribution 

in the boundary layer. 
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Arc length x 


Figure 35.- Velocity distribution of the potential 
flow for a wing profile. 



Figure 36.- Concerning the calculation of the friction layer on the 

circular cylinder. 


y 



Figure 37.- Wake flow behind the flat plate in longitudinal flow. 




Figure 39.- Asymptotic velocity distribution in the wake behind the 

flat plate in longitudinal flow. 





2 


1 


1 


2 


Figure 41.- 


? 


The velocity profile of the plane jet. 





Figure 41(a).- Boundary-layer profiles for the potential flow U(x) = u-jx 
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Figure 42.- Application of the momentum theorem for the flat plate in 

longitudinal flow. 



Figure 43.- Velocity distribution in the boundary layer on the flat plate 

in longitudinal flow. 

(a) Linear approximation. 

(b) Cubic approximation for the velocity profile. 



Figure 44.- Application of the momentum theorem to the boundary layer 

with pressure gradient. 
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Figure 45.- The universal functions F(y/6p) and G(y/8p) for the 
velocity distribution in the boundary layer according to Pohlhausen. 



Figure 46.- The one-parameter family of velocity profiles according 

to Pohlhausen. 
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Figure 47(a).- Auxiliary functions of the boundary layer calculation 
according to Holstein (cf. table 5); x and F( *) against k . 
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Figure 48.- Integration of the differential equation of the boundary layer 
according to Pohlhausen and Holstein (profile J 015; c a = 0). 


Joukowsky profile d/t = 0.15 


Figure 49.- Result of the boundary-layer calculation for the example 
according to figure 48 (profile J 015; c a = 0). 
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Stagnation point profile 

U m u f x 
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Figure 51.- Potential -theoretical velocity distribution on the elliptic 
cylinders with axis ratio a = 1, 2, 4, 8 for flow parallel to the 

major axis (A = laminar separation point), t' = half the circumference. 
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Figure 52.- Result of the boundary-layer calculation for the elliptic 
cylinders of axis ratio a^A^ = 1, 2, 4, 8. 








= Inflection point 
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Figure 54.- Potential -theoretical velocity distribution for the Joukowsky 
profile J 015 for c a = 0; 0.25; 0.50; 1.0. 











s/f-aoi9z 















Suction side 



Figure 57. - Potential flow and laminar friction layer for the Joukowsky 

profile J 015 at c a = 1. 










JoukowsKy profile 




Figure 59.- Position of the laminar separation point for the JoukowskM^v 
profiles of figure 58; suction side. 





e 60. - Position of the laminar separation point for the Joukowsky 
profiles of figure 58; pressure side. 
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Figure 61.- Pressure distribution and separation on a wing. 



Figure 62.- Boundary layer with laminar separation avoided. 


Ufxl 



Figure 63.- Potential flow with separation: U / < 0; U // < 0; 
sometimes without separation: U / < 0; U" > 0. 
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Figure 64.- Divergent channel. 



Figure 65.- Prevention of separation on wing by blowing. 



Figure 66.- Prevention of separation by a slotted wing. 



Figure 67.- NACA cowling for prevention of separation. 
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Figure 68.- Prevention of separation on wing by suction. 



Figure 69.- Flat plate in longitudinal flow with suction. 



Figure 70.- Asymptotic velocity profile on flat plate in longitudinal 

flow with suction (I) s* = v/-v 0 (II) 5 * = 1.73 -\f^- • 

/ O 0 


